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We have studied the Knight shift K(f, T) and magnetic susceptibility Jt(T) of heavy electron materials, modeled by the 
infinite (7 Anderson model with the NCA method. A systematic study of K(f, T) and x(^) f° r different Kondo temperatures 
To (which depends on the hybridization width T) shows a low temperature anomaly (nonlinear relation between K and x) 
which increases as the Kondo temperature To and distance r increase . We carried out an incoherent lattice sum by adding 
the K(r) of a few hundred shells of rare earth atoms around a nucleus and compare the numerically calculated results with 
the experimental results. For CeSn3, which is a concentrated heavy electron material, both the 119 Sn NMR Knight shift and 
positive muon Knight shift are studied. Also, lattice coherence effects by conduction electron scattering at every rare earth site 
are included using the average-T matrix approximation. The calculated magnetic susceptibility and 119 Sn NMR Knight shift 
show excellent agreement with experimental results for both incoherent and coherent calculations. The positive muon Knight 
shifts are calculated for both possible positions of muon (center of the cubic unit cell and middle of Ce-Ce bond axis). Our 
numerical results show a low temperature anomaly for the muons of the correct magnitude but we can only find agreement 
with experiment if we take a weighted average of the two sites in a calculation with lattice coherence present. For YbCuAl, the 
measured 27 A1 NMR Knight shift shows an anomaly with opposite sign to CeSn3 compound. Our calculations agree very well 
with the experiments. For the proposed quadrupolar Kondo alloy Y0.8U0.2Pd3, our 89 Y NMR Knight shift calculation doesn't 
show the observed Knight shift anomaly. 

74.70.Vy, 74.65.+n, 74.70.Tx 



I. INTRODUCTION 

Many heavy electron materials show Knight shift 
anomalies, which are a deviation from a linear relation 
of the Knight shift K(T) to the magnetic susceptibility 
x{T) below the Kondo temperature To. The origin of 
the Knight shift anomalies has been a subject of great 
interest in the condensed matter .community over a pe- 
riod of nearly twenty five years ErO. If the impurities in 
metals have local magnetic moments, they display inter- 
esting properties comparing to metals with nonmagnetic 
impurities, such as a resistivity minimum and anomalies 
in specific heat and susceptibility. This Kondo effect is a 
consequence of interaction between the magnetic ion and 
conduction electron. The central physical concept is that 
the many body screening cloud surrounding a Kondo im- 
purity site should give rise to an anomalous temperature 
dependent Knight shift at nuclear sites due to the cou- 
pling of the localjinnriLent to the nuclear spin through the 
screening cloucffloml Such a "non-linear Knight shift 
anomaly" is to be distinguished from the non-linear sus- 
ceptibility related to the field dependence of \- Another 
way to describe this effect is to say that in the absence 
of an anomaly, the contribution K(r, T) from a local mo- 
ment at distance r from the nucleus can be written as 
/(r)x(T). This factorization does not hold if there is an 



anomaly (instead K(r, T) = /(r, T)x(T) due to the tem^ 
perature dependent polarization cloud). After HeegerEJ 
suggested that the anomalous spin cloud was detected at 
low temperatures, the central question has been whether 
a conduction electron spin cloud with huge coherence 
length !;k = hvp/ksTo, where vp is the Fermi veloc- 
ity and To is the Kondo temperature exists. The main 
motivation of this paper is to help clear up this conflict. 
The Knight shift calculations presented here are the first 
performed using a realistic impurity model. 

This paper is organized as follows. In Sec. |l|, the 
Kondo effect and the Knight shift anomaly in heavy elec- 
tron materials will be reviewed. First, general character- 
istics of the Kondo effect are discussed. We will then 
review the history of Knight shift anomalies in heavy 
electron systems. In Sec. [II, the model Hamiltonian is 
introduced for both Ce and U compound. Also we review 
the methods we have used to evaluate the Knight shift 
(the Non Crossing Approximation(NCA) and average-T 
matrix approximation (ATA)). Our formalisms for nu- 
merical calculations are explained, and a detailed deriva- 
tion of the Knight shift Feynman diagram will be given in 
the appendices. In the next section, the numerical results 
for Ce and Yb ions which are single channel Kondo ma- 
terials and for U ions in a proposed quadrupolar Kondo 
alloy will be examined and compared with the exper- 
imental results. The calculated NMR Knight shift of 
Ce and Yb compounds show low temperature anoma- 
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lies and agree well with the experimental results. But, 
there is no calculated Knight shift anomaly for the pro- 
posed quadrupolar Kondo U compound, in contrast to 
experiment. The last section includes conclusions and 
directions for future work. 

II. REVIEW 
A. Kondo Effect 

The existence of localized moments in dilute alloys 
that couple to conduction electrons has important conse- 
quences for— the electrical properties. It has been known 
since 1930til that the resistivity has a rather shallow 
minimum occurring at a low temperature that depends 
weakly on the concentration of magnetic impurities in- 
stead of dropping monotonically with decreasing temper- 
ature like metals with nonmagnetic impurities. In 1963, 
Kondo£3 explained that this minimum arises from some 
unexpected features of the scattering of conduction elec- 
trons off a local magnetic moment, with the simplified 
model Hamiltonian 

ka kk' 

where cg Q is the annihilation operator of the conduction 
electron, a is the spin index, a are Pauli matrices, Si 
is the spin operator of the impurity, and J < is the 
exchange coupling. Kondo discovered that the magnetic 
scattering cross section is divergent in perturbation the- 
ory. The anomalously high scattering probability of mag- 
netic ions at low temperatures is a consequence of the dy- 
namic nature of the scattering induced by the exchange 
coupling and the sharpness of the Fermi surface at low 
temperatures. Subsequent analysis by Kondo and others 
has shown that a nonperturbative treatment removes the 
divergence, yielding instead a term in the impurity con- 
tribution to the resistivity that increases with decreasing 
temperature. In spite of the simple model Hamiltonian, 
a magnetic Si — 1/2 local moment interacting with the 
conduction electron gas, this result is an indication that 
the problem is explicitly a many body problem, meaning 
that the electron in state k which is being scattered is 
sensitive to the occupation of all other electron states q. 

For this single channel Kondo model, there is only one 
characteristic energy scale, the Kondo temperature To, 
provided that the temperature T is much smaller than 
the conduction electron bandwidth D, and corrections 
of order T/D are neglected. The Kondo temperature is 
given by 

k B T = D[N(0)J]^ 2 exp{-l/N{Q)J), (2) 

where iV(0) is the conduction electron density of states 
at the Fermi level. Any physical quantities are universal 



functions of T/Tq. At low temperature, with all material 
properties buried in To. 

This Kondo model can explain the anomalies in the 
transport coefficients, specific heat and magnetic suscep- 
tibility for some alloys with magnetic impurities. The 
Kondo effect is characterized by the development of the 
Kondo resonance peak with width of order of To. At high 
temperatures T > To, the impurity resistivity increases 
logarithmically as the temperature decreases, and satu- 
rates to a finite value at low temperatures below T < Tq. 
The magnetic susceptibility has a Curie- Weiss form at 
high temperatures and shows Pauli paramagnetism at 
low temperatures with %(T = 0) ~ 1 /To, concomitantly 
C/T(T = 0) ~ 1/T . The X behavior is explained by 
the fact that the magnetic moments which exist at high 
temperatures are screened out by the conduction electron 
spin clouds at low temperatures with the formation of a 
singlet ground state. This conduction spin cloud and the 
Knight shift anomaly will be discussed more in the next 
subsection. 

B. Review of the Knight Shift Anomalies of Heavy 
Electron Materials 

There have been many theoretical and experimental 
works about the Knight shift anomaly for the heavy elec- 
tron materials. Whether there is an observable conduc- 
tion electron spin cloud with huge coherence length has 
been an issue in condensed matter physics more than 25 
years and this is our main motivation to carry out this 
paper. 

In the simplest approximation, the added electron is 
bound into a singlet with the impuritytJ. Because of 
the falloff the amplitude of wave function with the en- 
ergy efe, only states within roughly 5k ~ (fesT)/ Ep)kp 
of the Fermi surface are involved in the singlet. As a 
result, in coordinate space, the singlet wave function ex- 
tends to a very large distance of order (<5fc) -1 . Heeger et 
al.t-H calculated the susceptibility at T = OK using the 
Appelbaum-Kondo theoryiilEj and found 

X = XPauli + XL + XQ (3) 

where 

XL = XQ = \JW 2 /(9/S)k B T } (4) 

Equations (^) and (||) give the very interesting result that 
one half the excess susceptibility is localized on the impu- 
rity site (xl) and one- half is associated with the partially 
polarized quasi particle (xq)- The associated spin polar- 
ization around the partially magnetized impurity is given 
by 

er(r) = <r + <7 RKKY (r) + cJ Q (r) (5) 
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where oo is the uniform polarization due tp-the external 
field, <Jrkky{t) is the usual RKKY termO~0 which r 
dependence is given by 



(l/r 3 )cos(2fc F r), 
and ctq (0) is the quasiparticle term with 



(6) 



< <W) >=< S z > ^ N £ 



3 AT k B To , sin k F r . 2 2 

) ln (>7f)- ( 7 ) 
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where = (2Ep /ksT^jk'p 1 . This expression is valid 
for r < and at greater distances erg(r) rapidly ap- 
proaches zero. In both the RKKY term and (Jq{t) term 
the value < Sz > is not the free spin value but is de- 
termined by the local susceptibility \l- The existence of 
the RKKY term for T <T was shown by SuhlE3. 

The technique of nuclear magnetic resonance has been 
of primary importance in the development of our current 
understanding of the localized moment problem. The 
reasons are twofold. First, the nuclei in the host metal in 
the vicinity of the impurity are sensitive to local pertur- 
bations in the spin density(via the hyperfme interaction) 
and the charge density(via the nuclear quadrupolar inter- 
action). Moreover, the nuclei themselves are only weakly 
coupled to the electronic system and therefore act as pas- 
sive "spies" into the phenomena of interest. Secondly, 
the nuclear relaxation is sensitive to the low-lying ex- 
citations of the electronic system and consequently can 
provide information on the dynamical aspects of the im- 
purity problem. For the most part NMR experiments in 
heavy electron compounds are carried out on the nuclei of 
the non-/ ions, so that coupling to the / moments occurs 
via indirect interactions such as transferred hyperfme and 
dipolar fields. „ 

The NMR experiment of Boyce and Slichtercl on Fe 
impurities in Cu metals showed no evidence for a Knight 
shift anomaly at low temperature and was interpreted to 
indicate the absence of this screening cloud or at least a 
screening cloud of size of the order of a lattice spacing. 
These mixed results have led to theoretical discussions 
about the size of any screeningpCflpduction electron spin 
cloud or even whether it existQtfflEj£3. 

In contrast, pronounced Knight shift anomalies have 
been observed in the concentrated heavy electron materi- 
als CeSi^M and YbCuAlErOl, which have been described 
as Kondo lattice systems with To > 400K. In view of 
the Boyce-Slichter result, the question is raised whether 
these anomalies represent a coherent effect of the peri- 
odic lattice rather than a single ion effect. However, re- 
cent experiments on the proposed quadrupolar Kondo 
alloyE£lEZI Yi_ x U x Pd3 demonstrate that for concentra- 
tions of 0.1-0.2 there are pronounced non-linearities in 
the Y Knight shift for sufficiently large distances away 
from the U ionEIl 

Ishiia calculated the field induced spin polarization for 
the degenerated Anderson model and confirmed that an 



anomalous spin cloud is formed outside of the Kondo 
screening length £r- at T — > 0. The spin polarization <r(r) 
for electron-hole symmetry case (e<j = — E7/2) is given by 



XH cos2fc F r 
cr(r) = ; — - — II 



9VB 



47J-J-3 



2E F 
k F rA 11 ' 



(8) 



where 7 = xtt for [/ ^ and 3/2(2/ + 3)xtt f° r U — » 00, 
where / is the orbital angular momentum. Xtl changes 
from ( ffA i B ) 2 /27rr for U = to (gfj, B ) 2 S(S + l)/3k B T for 
s — d Exchange model (U — 00). Therefore the coherence 
length £k = 2E F /k F Aj varies from 2hv F /T at U = 
to 2nhv F /3k B To in the s — d limit. Also, for the strong 
U /ttT limit, spin polarization is calculated for r < It 
is given by 



a(r) 



4r xH cos2k F r 



(9) 



irU g[i B 47rr 3 

By the relation W 2 /U = - J/2JVB, the above spin po- 
larization is just the RKKY contribution. Comparing the 
equations (|j) and (||) , the conduction electron spin cloud 
which is formed outside the Kondo coherence length has 
the RKKY form but with 2N/N(0)J times bigger am- 
plitude than the spin cloud inside the Kondo coherence 
length. 

Chen et alB calculated the zero frequency response 
function C (r, T) around a magnetic impurity, using a 
perturbative thermodynamic scaling procedure and non- 
perturbative renomalization group method for S = 1/2 
Kondo model. The host nuclei near magnetic impurities 
at positions r displays satellite resonances in the tail of 
the main magnetic resonance signal, with a Knight shift 
given by K + AK(r), where K = pog e /2 is the Knight 
shift of the pure host. Then AK is 



AK = C(r, T)-K 



(10) 



Chen et alu showed that this Knight shift is factorized 
into a product of temperature and spatial dependent 
functions, specifically 



AK = F(r) X (T) 



(11) 



where x(T) is the magnetic susceptibility. The computed 
F(r) has the RKKY form. 

The conduction electron spin polarization in the An- 
derson model has been studied with the NCA previously 
by W. Pollwein et alu. However, this study was carried 
out only for the spin 1/2 model with infinite Coulomb 
repulsion, and for a limited parameter regime (only re- 
sults for very low To values and short distances (r <!C S,k) 
have been numerically calculated). In consequence, no 
strong evidence was found for a Knight shift anomaly in 
this previous work. _ _ 

Recently, E. S0rensen and I. AffleckcJLJ showed that 
the Kondo coherence length £ F = hv F /ksTg, varies 
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when temperature changes by combining a finite size scal- 
ing ansatz with density matrix renomalization group cal- 
culations. They write the scaling hypothesis for three 
dimensional susceptibility is 

XimpM = * - P/2 

cos(2fc F r) 
= -^2^f(rT/v F ,T/T ) 

(forr> l/k F ;T,T ^:E F ), (12) 

where / is a real universal scaling function, p/2 is the 
standard Pauli bulk susceptibility with p the density of 
states per spin. At higher temperatures T > Tq, the 
local susceptibility shows RKKY behavior and at lower 
temperatures T < Tq it has a local Fermi liquid form. 
So the Knight shift has longer range at low temperatures 
where the conduction electron screening cloud has formed 
than at high temperatures where it has not. Sorensen 
and Affleck Jried to explain the experiment by Boycc 
and SlichterB by the possible factorization of the scal- 
ing functions deep inside the screening cloud where the 
experiment was done. 

III. MODEL AND FORMALISM 
A. Model Hamiltonian 

In our work, we use the on-site Coulomh interaction 
U = oo single impurity Anderson modelEj. The An- 
derson model can be canonically transformed by the 
Schrieffer-WolfE-transformation to the Kondo model at 
U = oo limited, and is a good model Hamiltonian 
to describe heavy electron materials with / electrons. 
Schematics of this model are shown in Fig |l| for Ce ions. 
This U — oo model can be a good approximation in the 
limit when the ratio of the virtual 4/ level width T to the 
on site Coulomb repulsion U is small. For real materials 
this interaction energy is of the order of U — 5eVt3'E2l, 
and the hybridization width is of the order of T = O.leV. 
For Ce(Yb) ions, we keep /° and / 1 (/ 14 and / 13 ) config- 
urations, and for U ions, which are proposed two channel 
quadrupolar Kondo alloys, we keep f 2 and / 3 configura- 
tions. For Ce and Yb ions, spin orbit coupling is included 
and for U ions, spin orbit coupling and also crystal field 
splittings are included. The crystal electric field(CEF) 
will split the spin-orbit multiplet and can mix two dif- 
ferent angular momentum multiplets(j, j'). But the cor- 
rection of the mixing term between two different j's by 
CEF is small, we only consider the splitting effect. These 
crystal electric field effects in the Anderson model were 
first considered by Hirsto based on group theory and it 
will be discussed more in the Appendix [b|. 

We shall first discuss the situation for Ce 3+ and Yb 3+ 
ions, and write down the model only for the Ce case(the 



Yb 3+ ion has a lone 4/ hole and our procedure describes 
this with a simple particle hole transformation). 
For a single Ce site at the origin, the model is 

H = H C + Hf +H cf +H Z (13) 

with 

ka 

the conduction band term for electrons with a broad 
featureless density of states of width D, taken to be 
Lorentzian here for convenience, with 

Hf = '£e fj \f 1 jm j ){f 1 jm j \ (15) 

where j = 5/2, 7/2 indexes the angular momentum mul- 
tiplets of the Ce ion having azimuthal quantum numbers 
to, with e /5/2 = -2eV,e f7 / 2 = e/5/2 + A so = -1.71eV 
(we take the / configuration at zero energy), with 

«c/ = E I Vi^/jfOHf^miil + h.c. } (16) 

kjrna 

where V^ jam . = VY imj _ a (k){Zmj - a, l/2a\jm j )/\/W a , 
V being the one particle hybridization strength and N s 
the number of sites. We can rewrite the hybridization 
Hamiltonian as 

Wc/=£[ V k c\. mi l/OQX/V^I + h.c. ] (17) 

kjrrij 

where 

C lm 3 = J2^ m l a ^ m ^ C lma ( 18 ) 
mer 

with Ckma — k J dkY£ m {k)cz . For the Yb case, the hy- 
bridization Hamiltonian is given by 

Wc/=E[ VkC hrn 3 \f U J - ^X/ 14 0| + h.C.}. (19) 

kjrrij 

H z , the Zeeman energy of the electronic system for a 
magnetic field H z applied along the z-axis is given by 

U z = -usH^^an^ -^9jrn\f l jm){f l jm\ ]. 

ka i' m 

(20) 

In addition to this, we must add a term coupling the 
nuclear spin system to the conduction electrons, which 
we take to be of a simple contact form ~ 1(f) ■ S(f) 
for each nuclear spin 1(f) at position f with S(r) the 
conduction spin density at the nuclear site, and a nuclear 
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Zeeman term. In terms of the parameters, the Kondo 
scale characterizing the low energy physics is given by 



k B T = D{ 



r )1/N a( D 



\N B */Ng 



exp( 



N a T> 



(21) 



where the single particle hybridization width T — 
irN(0)V 2 with N(0) which is the density of the states 
at the Fermi energy. Other parameters are defined in 
Table |. 

For the Yx- x U x Pd^ which has the cubic AuCu3 struc- 
ture, the crystal field effect (CEF) must be included. 
This crystal electric field effects lift the angular momen- 
tum degeneracy of U ions and their spin-orbit multiplet 
decomposes into irreducible representation of the cubic 
field. The f 2 , j = 4 Hund's rule ground state of U com- 
pound is split to a T3 nonmagnetic doublet, and T4 
magnetic triplet and Ti singlet statesu. And / 3 , j = 9/2 
spin-orbit multiplet is split to T 7 doublet and two r 8 
quartets. In our calculation we choose j = 4 T 3 for the 
ground state for f 2 configuration and j — 9/2 for the 
ground state of the f 3 configuration. Fig. || shows the 
schematic configuration diagram. All parameter values 
are listed in Table || in the unit of D. For an explicit 
derivation of Hamiltonian for U ions, see the Appendix^. 



B. Non Crossing Approximation 

We treat the Anderson Hamiltonian with the non- 
crossing approximation (NCA), a self- consistent dia- 
grammatic perturbation theory discussed at length in 
the paper of Bickers et alx3. This is useful because 
this method provides ways of calculating the dynamic 
response functions, such as the one electron Green's func- 
tions and dynamic susceptibility and it makes possible a 
more extensive comparison between the theoretical pre- 
dictions and experimental results. 

In the NCA, we do the 1/N expansions with the new 
variable N, the large orbital degeneracy of the ground 
state of / electrons. These simple approximation schemes 
work very well for values of N that are of interest in 
applications to rare earth impurities. For example, the 
lowest spin-orbit split multiplet for Ce 4/ 1 has j = 5/2, 
corresponding to N = 6, and for Yb 4/ 13 , j = 7/2, 
corresponding to N — 8. Even for N — 2 one can get 
good semi-quantitative results, which can be accurate to 
within a few percent for some quantities. 

It was noted that, in a perturbative 1/N expansion, 
the /-electron spectral density exhibits a singularity at 
lu = To, with To ~ _Dexp(7re//AT). This singularity re- 
mains order by order, preventing a complete description 
of / photoemission and electronic transport. In order 
to remove this singularity, it is necjessariz-to perform an 
infinite-order resummation in 1 /N Ea^jTcZl. 



In the NCA , our starting basis is the conduction band 
plus the atomic Hamiltonian projected to the atomic elec- 
tron Fock space and treat the hybridization between the 
conduction band and the atomic orbital as a perturba- 
tion. The strength of this approach is that the strong 
on-site Coulomb interaction for atomic electrons is in- 
cluded at the outset. The conventional Feynman diagram 
technique which uses Wick's theorem can not be applied 
for strongly correlated problems with restricted Hilbcrt 
spaces. Pseudo particle Green's functions are introduced 
for each atomic electron occupation state which is nei- 
ther fermionic nor bosonic (i.e., f 1 j — 5/2, 7/2 and /° 
in the present model for Ce ions). The pseudo fermion 
Green's functions for f 1 j — 5/2, 7/2 angular momentum 
multiplets are 



Gjm(z) 



1 



Z 6 fj ^jmi 2 ) 

and the pseudo boson Green's function for the / is 
D(z) = ' 



z - S (z) 



(22) 



(23) 



Then we insert a self-energy into the propagators of 
pseudo particles. This gives coupled integral equations 
for the ionic propagator self energies, So(z), T,j m (z). 
^From the leading order diagrams of Fig. ||, the coupled 
equations for the self energies are 



£o(z) = ^ 2 $>^ 



fk 



■f z + e k - efjm - T, jm (z + e k ) 



X jm (z) = V 2 J2 



l-fk 



z-e k - E (z - £fc) 



= V 2 J2(l-fk)D(z-e k ), 



(24) 



(25) 



where V is the hybridization strength between the con- 
duction band and the atomic orbitals and Nj is the de- 
generacy of the spin-orbit multiplet j. It is convenient 
to introduce the spectral functions A,- m (u;),and B(uj) for 
pseudo-particle Green's functions. 



Aj m {w) 
D(iv) 
B(v) 



dp Aj m (p) 



-Im Gj m (uS) 
dC B(Q 

7T %V — £ 

Im D{v) 



(26) 



(27) 



In addition to spectral functions Aj m {p) and B{C,), it is 
necessary to introduce negative frequency spectral func- 
tions Aj^(p) and These spectra are given by 



5 



A^(p,T) = A jm (p,T)e-^- E ^ 



(28) 
(29) 



where Eq is the ground state energy relative to the Fermi 
energy. The impurity partition function Z^f is given by 

jm 

= e -PE J ^ { J2a^(c,t)+b^(c,t)}. (30) 



At T -> 0, Z 4 /(T) becomes 



Z 4/ (0) 



(31) 



The iteration of these coupled equations for the self 
energies generates a set of diagrams which includes all 
non-crossing diagrams, but does not correspond to any 
specific order in the 1/N expansion by treating NV 2 as 
0(1)) where V is the hybridization strength between the 
conduction electron and the atomic orbitals. The set 
of diagrams summed by these equations includes all the 
terms of order 0(1) and 0(1/N) and a subset of con- 
tributions from the higher order terms. The lowest or- 
der skeleton diagrams which are not included are of or- 
der 0(1/N 2 ). All the diagrams that enter at 0(1) and 
0(l/N) have non-crossing conduction lines. Specifically, 
the leading order vertex corrections, which are 0(1/N 2 ), 
are not included in the NCA. 

These self-consistent integral equations are solved to 
second order in the hybridization for the ionic propaga- 
tor self energies. Then physical properties, such as the re- 
sistivity and magnetic susceptibility x> are calculated as 
convolutions of these propagators. Fig. [| shows a leading 
order Feynman diagram for the static magnetic suscep- 
tibility and its convolution integral is given by Eq. (32) 
This is discussed more in section 
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The NCA shows a pathological behavior (due to the 
truncation of the diagrammatic expansion) for a tem- 
perature scale T p <C To in this conventional Anderson 
model. However, provided the f 1 occupancy n/ > 
and N > 4, this is not a problem, as shown in Ref.E 
in that comparison of NCA results with exact thermody- 
namics from the Bethe-Ansatz shows agreement at the 
few percent level above T p . Hence, this is a reliable 
method for our purposes. 

Our numerical procedure, briefly, consists of solving 
the NCA integral equations for the Anderson Hamilto- 
nian specified above on a logarithmic mesh with order 600 
points chosen to be centered about the singular structures 
near the ground state energy E ~ e/s/2- We then feed 
the self-consistent propagators for the empty and singly 
occupied orbitals into the convolution integrals obtained 
from the diagram of Fig. ||, which allows for evalua- 
tion of the Knight shift at arbitrary angle and distance 



from the nuclear site. This will be explained more in sec- 
tion III D . It is convenient to take the nuclear site as the 
origin in this case leading to phase factors e~ lk ' R in the 
hybridization Hamiltonian 7i c f , where r is the nuclear- 
Ce site separation. These factors give the oscillations and 
position space angular dependence in the Knight shift K. 



C. Magnetic Susceptibility 

The static magnetic susceptibility is a direct indicator 
of the nature of the ground state for the Kondo and An- 
derson model. Near room temperature, the susceptibility 
x(T) displays a Curie- Weiss temperature dependence for 
most heavy-electron materials, x is linearly related to the 
Knight shift K(T). At low temperature, x(T) does not 
follow the Curie- Weiss law and for some systems, a linear 
relation of the Knight shift K(T) to the magnetic sus- 
ceptibility breaks down and hence shows a Knight shift 
anomaly. 

In general the magnetic susceptibility comes from self- 
correlation of the conduction band magnetization(x cc ), 
the self-correlation of the / magnetization^//), and the 
mutual correlation of / and band components(x c /). The 
leading diagram comes from the second term where the 
only / electrons are coupled to the field and this gives 
a good approximation to the overall susceptibility. This 
diagram for magnetic susceptibility is in Fig. ^. Then, 
in the NCA, the magnetic susceptibility in the zero field 
limit can be writtei£3 



X(T) = £ 



Z4f(T)' 



A\2((,T)BeG jm (t,T) (32) 



where fij is the effective magnetic moment which is de- 
fined as /j, 2 = (gjHB) 2 j(j + l)/3 where gj is the Lande 
g factor for the j multiplet and Z 4 /(T) is the impurity 
partition function. 

Also we can get the van Vleck magnetic susceptibility 
between j and j' angular momentum multiplet. 

" f \ _ 2 1 (jrnj \J Z + S Z \j' mj , ) I 2 8 



j ^ Ajmj (C + uj)A^.,(0 



(33) 



With (Z,m\J z + S z \l,m/) = — (V49 - Am 2 /U)5 mm >, the 
total van Vleck susceptibility for j = 5/2 and j' = 7/2 is 

nij ,rrij/ 

-r- [1 z£v I v A - (c+w)A ^ (c) - (34) 



G 



The susceptibility sum rule is derived from the zero fre- 
quency limit of the Hilbert transform of x vv (uj 1 T), 



(35) 



In our calculation, both the j — 5/2, 7/2 contribu- 
tions to the static magnetic susceptibility and the van 
Vleck susceptibility are considered for Ce ions. For Yb 
ions, only the static susceptibility of the j = 7/2 ground 
spin-orbit multiplct is calculated because that the en- 
ergy gap between two spin orbit multiplets (j — 7/2 and 
j = 5/2) is large ( about 1.2eV) For Yi_ x U x Pd 3) only 
the van Vleck susceptibility between different T states 
is calculated because the assumed T 3 ground state is a 
non-magnetic doublet. 



D. Knight Shift 

Knight shift measurements on the nuclear spins of non- 
/ ions in Kondo or heavy electron materials can probe the 
local induced magnetic fields. The additional fields come 
from all the possible polarization sources, such as con- 
duction electron spin polarization. For / electrons, the 
radius wave function is small and they are well screened, 
so there is little possibility of direct overlap interactions 
between the nuclear spins and local moments. In par- 
ticular, the polarization of conduction electrons by the 
polarized local Kondo impurities, i.e. the transferred 
/-electron polarization, is usually expected to have the 
most significant temperature dependent contribution. 

The Knight shift of Heavy electron materials is induced 
by the indirect interactions of the magnetic impurity and 
host nuclear spin mediated by the conduction electrons. 
Without the charge fluctuations introduced by the hy- 
bridization interaction between conduction electron and 
/ electron, this indirect interaction has the RKKY inter- 
action form. So at high temperatures, the Knight shift 
follows the RKKY interaction and at low temperatures 
where the Kondo effect appears, the Knight shift can 
show deviation from the RKKY form. 

To calculate the Knight shift we need to evaluate the 
Feynman diagram in Fig |5| which is the lowest order dia- 
gram coupling the nuclear spin to, say, Ce magnetic mo- 
ments, ignoring direct Coulomb exchange coupling. All 
the propagator symbols are explained in the figure. For 
the incoherent calculation, the conduction electrons are 
assumed to belong to a broad, featureless, and symmetric 
band of half width D. The conduction electron propa- 
gator is taken to be a bare electron propagator, i.e. it 
includes no self energy effects reflecting multiple scatter- 
ing off the /-sites. When the the lattice coherence effects 
of conduction electrons are included, the self energy aris- 
ing from the scattering of conduction electron at every / 



site is included in the conduction electron propagator in 
an approximate way. 

The Knight shift for a nuclear spin or muon at f is 
approximated as 

8nV 2 g e p 2 B 



K{r,9) 



£/«'(M) / 



3N S n 

x / k' 2 dk'I jf (e k ,e k >), 



k z dk 



(36) 



where g e is the electron g-factor and \ib is the Bohr mag- 
neton. Here 9 is the angle between z axis and the bond 
direction r which connects the nucleus or muon to a given 
/ ion in the crystal and a is the angle between the field 
axis and bond axis. Ijf (efc, ey) is given by 

Ijj'( e k, £fc') = j3(kr)j 3 (k'r)—— V* G (ka, iu))G (kW , kS) 
xGj m (i(u + v))Gji m i{i(u + v))D(iv) 

+(l-/(e fe ))I 2 ^(e fc )), (37) 

where 

2 [d( [dp (B^(QA jm (p)ReG fm ,(p) 
lijj'{ek) = / — / — 

Z if J 7T J 7T E/c+C-P 

= ^-/^B(-)(0RfiG im (C + 6 fc ) 



xReGfm'iC + efc) 



(38) 



hjj'{ek) 



2 [dp fdCB( QA ( j J(p)RcG f „Ap)) 
£k + C~ P 



Z/± f J 7T J TT 

J_ fdp A (-) 

The derivation is explained in detail in Appendix |c[ 
We can analytically evaluate the inner k' integral as 



A^(p)ReG ym ,( P )ReD(p-e k ). (39) 



J2 ,,j Mk'r) 



k dk' 



k 2 - k! 2 ' 



(40) 



and the results are presented in Appendix 

For a magnetic field in the z direction {i.e., 9 — a), 
and for / = 1/2, the angular dependent function in the 
Knight shift, fjj'(ct,8), is given by 

/*'(») = cos2 °& + ^(4" + fjr + /# + 

(41) 

Here f? f {6) is 

rrij ,Trij/ 

= £ a(jrrij\J z + S z \j'm jl )(j'm j/ \3m 3 ; -a) 



rrij ,0,1713 



x(3m 3 ; ■xa\jm j )\Y 3m3 {f)\ 2 . 



(42) 
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and and/2±(0) are 



rrij rrij/ a,/3 

XI {h , m j \J±\j',m r ){j'mj'\Zm'^ hat) 

rrij ,rrijf ,7713 ,mg,a,/3 

x(3rn 3 ; 2«ka,^|3m 3 ; ^/3)(3m 3 ; ^Pljrrij) 
xY 3 * m3 (r)Y 3 m> 3 (r) (43) 
fjP( d )= Um j \S±\j\m j '){j , m j >\r)a*p{f > \j,m j ) 

X 0'>» 7i ji^±b' / :' 7i J')0'' m i'l 3m 3;2 a ) 

x(3m 3 ; -a|crj >/3 |3m 3 ; i/3)(3m 3 ; ^|jmj) 
The explicit values for fji (9) are 



6 sin 4 1 



/||(0) = — [2 + 12sin 2 0-9sin 4 



3 sin 4 01 



(44) 

(45) 
(46) 
(47) 



To assess the relevance of this single site physics to 
the periodic compound CeSn3, we also carried out inco- 
herent lattice sums over a few hundred radial shells of 
Ce atoms around the Sn nucleus. For each atom in a 
shell the distance r and angle 8 of Ce ion is calculated 
and the Knight shift of each ion is evaluated for given 
position. Then the contribution from each ion is added 
to get the total Knight shift. For the Y^^U^T^, im- 
purity configuration averaging is also carried out. This 
single site physics is known to be a good approximation 
at high temperatures where the ions are incoherent with 
one other, and known to provide a very accurate descrip- 
tion of the thermodynamics in many cases. For CeSn3, 
given the tetragonal symmetry at the Sn site, we fixed 
the field in the z direction and averaged over the xz,xy 
and yz host planes for the Sn nucleus. Note that yz plane 
is equivalent to xz. 

For YbCuAl, which has hexagonal symmetry, we have 
to consider three possible field directions, along the x, y 
and z axes. Al has an I = 5/2 nuclear spin and the NMR 
shift was obtained from derivative spectra of the central 
(1/2 -M- -1/2) NMR transition. Then 



B 



Uf (a, 9)= A cos a/?-, (8) + - sin a[f+ (9) + f~, (9)} 



'33 



Then, 



A= - (25 cos 5 9 - 26 cos 3 9 + 5 cos 9) 
B = - sin 8(25 cos 4 0-14 cos 2 0+1). 



fn(a,6) = — [2Acosa(l + 3 sin 2 6 

2 2 Y7T 

+Bsina(5 + 3 cos 2 0)], 
/ss(a, 9) = -^—[Acosa(l - 4 sin 2 6 

2 2 287T 

-Bsina(3 - 2 sin 2 0)]. 



(49) 



(50) 
(51) 



For a detailed derivation, see Appendix ||. 

To the extent that the dynamics of the empty orbital 
can be neglected, the Knight shift expression (Eq. (pr^)) 
factorizes into a nearly temperature independent RKKY 
interaction (modified due to the spin-orbit coupling and 
anisotropic hybridization from the original form) times 
the /-electron susceptibility. Thus, no anomaly results 
from the diagram in this limit. In this limit, the sus- 
ceptibility in the diagram corresponds to the the lead- 
ing order estimate used in RefEj to compare with exact 
Bethe-Ansatz results. 

We can gauge the effects of charge fluctuations with a 
simple approximation eS'EI. For T = 0, the empty orbital 
propagator may be written in an approximate two-pole 
form, one with amplitude 1 — Z, Z = irks^/NT, cen- 
tered near zero energy, and one with amplitude Z cen- 
tered at e g ~ et — kgT which reflects the anomalous 
ground state mixing due to the Kondo effect. The singly 
occupied propagator has a simple pole structure. 



B{uj) = (1 - Z)5(lu) + ZS(lu - e g ) 
A m (uj) = S(oj - ef) 
B(~\u) =Z5{uj-e g ) 



(52) 



Then 



h{e k ) = Z/(k B T - e k 
1 



Z . 



. Z.l-Z 
h(e k ) = -jr^ri + - 

KB-LO £k-Cf Cfc 



(53) 



Now, we can perform the k integral putting above equa- 
tions in Eq. (|37|). At low temperatures only conduction 
electrons which have momentum close to kp participate 
in the interaction. We can rewrite the radial momentum 
k as 



k = kp + (k — kp) fa kp + k 

nvp 



(54) 



(48) 



Then we can write 



where ff-, = fjp + ftp . Also A and B are given by 



r2± 



e ik F r e ie k r/hv F 



(55) 
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^From this we see that small €k gives large contribution 
to the Knight shift. The contribution from the integral I\ 
depends on whether r 3> hvf /UbTq. For r ^> hyf/ksTo 
only very small tk < fcsTb contribute the k integral and 
we can approximate l/(fcgTb — e^) 1/fcgTo and this 
contribution has the amplitude Z/ksTo = tt/NT. For 
r <C hvf //ebTo, the Knight shift has contributions from 
£fc > fcaTo and this term has the amplitude Z/D. The 
amplitude of the Knight shift outside of the coherence 
length £fc = TiVF/kBTo is D/ksTo times bigger that that 
at inside the the Kondo screening cloud. This term can 
give the Knight shift anomaly. Because e/ is much big- 
ger than the the first term of I2 gives conventional 
RKKY oscillations modulo the anisotropy and altered 
range dependence induced by the m, k dependence of the 
hybridization. The amplitude of the second term goes to 
zero above the Kondo temperature. This term also may 
contribute to the anticipated anomalous Knight shift, 
and within such a two-pole approximation may be seen to 
be finite within £k, have a stronger distance dependence 
in that regime, but possess an amplitude of order Z only 
within this distance regime. Beyond the amplitude 
is of order 1 /N and the shape of the spin oscillations is 
the same as that found from the high frequency pole of 
the empty orbital propagator. 

S0rensen and AffieckcjEa have noted that for a single 
impurity an additional contribution is present below the 
Kondo screening cloud which is not present in this calcu- 
lation. This corresponds to the diagrams shown in Fig. |^, 
in which scattering occurs off of the impurity for one or 
other conduction legs. This process will occur for any im- 
purity in a metal, and at a fundamental level corresponds, 
to the contribution induced by the field dependence of 
kF for the two different spin branches. The contribution 
then goes as 1/R 2 , but only outside the Kondo coher- 
ence length where the low temperature screening cloud 
can be regarded as a potential scatterer. (The 1/R 2 fol- 
lows trivially from differentiating the Friedel oscillation 
with respect to kF-) This contribution is of potential im- 
portance for any dilute system, but we argue that it is 
not important in our lattice context (or, for that mat- 
ter, for any system with an appreciable concentration of 
impurities). The reason is that the T- matrix insertions 
of Fig. H will go over to self-energy insertions in the lat- 
tice, as we sum over all possible sites. These self-energy 
insertions will simply provide the renormalized Pauli sus- 
ceptibility contribution to the Knight shift, which is not 
the dominant contribution. 



E. Coherent Lattice Effects 

Coherent lattice effects are included within the local 
approximation ( d = 00 limit) to the lattice model. In 
this approximation a conduction electron self energy is 



included using the average -T matrix approximationE3 
which assumes that conduction electron scatters off ev- 
ery /-electron site. This corresponds to a first iteration of 
the local approximation. In contrast, these multiple scat- 
tering processes are ignored in the incoherent limit. The 
NCA approximation treats intra-site interactions to all 
orders. In this calculation we consider intersite coupling 
which involves simple hopping process in perturbation 
theory (ATA) and ignore intersite interactions which in- 
volving transfer of particle-hole pairs between sites. This 
coherent lattice effect may reduce the Kondo screeinig 
lengthen. 

The Anderson lattice model for spin L/2 has a conduc- 
tion electron Green's function given byE3 



G c (k, to, T) = (w - e fc - E c (w, T))~ 1 

V 2 

= {0J- e k - 



)~\ (56) 



Where £™*(fc, u>,T) is the /-electron self energy arising 
from / — / interactions. The same results follow for the 
N-fold degenerate model. We remove wave- vector depen- 
dence of the self energy by neglecting the intersite inter- 
actions. The / electron self energy comes from the / — / 
interaction (£™* ) and hybridization with the conduction 
electrons. This hybridization energy is given by 



V 2 

V 1 N K *-< 



1 



(57) 



For a featureless symmetric band D(oj)V 2 becomes, for 
oj — > 0, ir, where T — 7riV(0)T^ 2 is the single particle hy- 
bridization width(7V(0) the density of state at the Fermi 
energy). Here the conduction electron can not be scat- 
tered at the same / electron site twice and this site re- 
striction gives the cancellation of the hybridization self 
energy term (— iT) in the full 4/ Green's functiorH. 

Specifically the band electimL self energy is written 
within this approximation as E3o 

E c (u,, T) = V 2 [{G if {u, T))(- x ) + V 2 D(w)](~ r> 

= V 2 G { :f(u,T), (58) 

where G4 / is the the full 4/ on-site Green's function given 
byS 

Ga/jM = 4~ I —[B(-Ht)G jm {cj + Z + i6) 

Zif J 7T 



Z i} 



(59) 



within the NCA. 

Since vertex corrections are neglected in the NCA, the 
conduction electron Green's function is determined com- 
pletely through Sj by hybridization and any resistance 
solely arises from the damping of band states due to the 
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imaginary part of / electron self energy. A realistic esti- 
mate of S/ is important to study lattice coherence effects. 
By the standard Fermi-liquid phase space arguments, 
the imaginary part of the exact, on-site f-electron self 
energy is given , for low frequencies and temperatures, 
by 



(65) 



ImX^/ (oj, T) 



(7Tk B T) 2 +T 



(60) 



In the NCA calculation, due to the approximation in- 
volved, the minimum value in — ImX^y (u>, T) does not 
occur precisely at the Fermi energy (though it differs only 
by a small fraction of T ) and is not equal to Itil. So in 
our numerical calculations, Im£4y(w,T) is extrapolated 
to u>, T — > and — Im£4/(0, 0) is replaced by T. 

In our calculation, the /-electrons have j — 5/2 and 
j = 7/2 states by spin-orbit coupling. Thus the conduc- 
tion electron self energy has two terms from each j state, 
viz. 

£ c ( W ,a,T) = IO- ! m i |^(fc)|fe,a)| 2 Gi7; ) (a;) 
= 4nV 2 J2\(j^j\3m 3 ^a)\ 2 



x\Y 3m3 (k)\ 2 G [ :;*\u) 

V l d(jr 4/5/2 + 4(jr 4/7/2i- 



(61) 



where the one particle hybridization strength V(K) = V 
is, taken to be independent of |fc| in this calculation, and 
we used 



Y,\(lm\3m 3 ^a)\ 2 \Y 3m3 (k)\ 

m 

Y^\(l, m \3m 3 ^a)\ 2 \Y 3m3 (k)\ 2 = - 



3 

4tt 



(62) 



With the inclusion of lattice coherence effects, the 
Ijj>(ek,£k') term in the Knight shift calculation is 
changed from the incoherent form, Eq. (p7j), to 



k 2 dk I k' 2 dk'I ir {e k ,e k r) 



ImJ(0ReJ(C)[/(0^(0 + (1 - f(0)h 3 f(0} 

(63) 



where 



j(0 = / dk 



k 2 j 3 {kr) 



(64) 



2 fdC fdp(B(-\C)A jm 0>)^Gj'm'{p) 





' 7T 






Z^fJ 


7T 



^4/ J 7T 



^- / ^0)ReG,, m , ( P )RcD(p - 0. (66) 



IV. RESULTS 

In this chapter, the numerically calculated results for 
the Knight shift K(T) and magnetic susceptibility %(T) 
of heavy electron materials such as CeSn3, YbCuAl, and 
U0.2Y0.8Pd3 will be presented. CeSn 3 and YbCuAl are 
concentrated heavy electron materials and U0.2Y0.sPd3 
is a proposed two channel quadrupolar Kondo heavy elec- 
tron alloy. 

First, the Knight shifts K(r,T) are systematically 
calculated for different values of the Kondo tempera- 
ture To which is controlled by the hybridization width 
r = ttN(0)V 2 , where N(Q) is the conduction electron 
density of states at the Fermi energy and V is the one par- 
ticle hybridization strength, for Ce and Yb compounds. 
These results show that the magnitude of the Knight 
shift anomaly depends upon the distance r between the 
local magnetic moment and the nucleus and the Kondo 
temperature To. There is an anomaly even for the small 
distance r. The magnitude of deviation between a lin- 
ear K .vs x relation is systematically increased when the 
distances are increased and the Kondo temperatures are 
increased. These results are shown in Fig. [ll] and Fig. |lj. 
These calculations support Ishii's idea, of an anomalous 
conduction electron spin density clouds which sets in be- 
yond the Kondo screening length = %vp /ksTo, where 
vp is the Fermi velocity. 

The lattice sum is carried out over a few hundred shells. 
In subsection IV E| the results for CeSn3 are compared 
with experiments. Both the 119 Sn NMR Knight shift and 
/isr Knight shift are studied. Also the influence of lattice 
coherence of the conduction electrons on both the NMR 
and positive muon Knight shift is investigated using the 
average T-matrix approximation and the numerically cal- 
culated results for CeSn 3 are compared with the exper- 
iments and also the calculated incoherent Knight shift. 
The calculated 119 Sn NMR Knight shift agrees well with 
the experiment. The incoherent psv Knight shift shows 
an anomaly but has opposite sign. For coherent case, the 
Knight shift from different muon site gives an anomaly 
with opposite sign. We may fit the experimental results 
by averaging out two psr Knight shifts. 

For YbCuAl, because of its complicated crystal struc- 
ture, the incoherent lattice sum is carried out over several 
thousand atoms. The calculated 27 A1 NM R Knight shift 
results are mentioned in subsection IV F. These results 



show excellent agreement with the experiments. 
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In the last subsection, the Knight shift and the mag- 
netic susceptibility of U0.2Y0.sPd3 is discussed. We do 
a full incoherent lattice sum and impurity configuration 
averaging for U0.2Y0.sPd3. Our study doesn't show the 
low temperature Knight shift anomaly like experiment. 



A. Systematic Calculations 

To see the systematic behavior of the magnetic sus- 
ceptibility and the Knight shift, we have calculated x(T) 
and K(r,T) for different Kondo temperatures To. For 
spin-orbit coupling and zero crystal field splitting, To is 
given by 



k B T = D(- 



*"l e /l 



■N g r' 



(67) 



where N g is the degeneracy of the ground spin orbit mul- 
tiplet, N ex is the degeneracy of the excited multiplet, tf 
is the energy level position of the ground multiplet, A so 
is the energy gap between two spin orbit multiplets, and 
r = ttN(0)V 2 is the hybridization width (D is the phys- 
ical Lorentzian bandwidth of the conduction electron). 
The values of the parameters we used for Ce and Yb ions 
are listed in Table |. 

For several different hybridization widths T, the 
Knight shift was studied as a function of temperatures 
and distance r between the local impurity spin and nu- 
cleus at fixed angle 9 between f and quantization axis 
k — k' . In these calculations all other variables such as 
e f , and A so were fixed to the values which give the best 
magnetic susceptibility x(T) fit to experimental results 
of CeSn 3 for Ce ions and YbCuAl for Yb ions (see sub- 
section IV B and subsection IV F for the parameter val- 



ues.). The Knight shift is scaled to the susceptibility by 
matching at high temperatures. For all the calculations 
the conduction electron band width D was assumed to 
be 3eV. Because of the small gap between j — 5/2 and 
j = 7/2 states of Ce compound (A so = 0.29eU for Ce 
ions and A so = 1.30eU for Yb ions) , the van Vleck term 
is included for only Ce compound studies. 

Fig. H shows the calculated Knight shift K{r, T) as 
a function of separation r and temperatures T at fixed 
angle 9 = and hybridization width T = 0.050667D 
(To = 430^) for Ce ions. We use a dimensionless vari- 
able kpr with the Fermi wave vector kp — 0.65A -1 in- 
stead of r. The Knight shift on a fine scale is shown 
in Fig. 0. The Knight shift shows an oscillatory RKKY- 
like behavior and the the total magnitude decreases when 
distance is increased. |A^(T)| at fixed distance r is first in- 
creased and then decreased (it becomes almost constant), 
when the temperature is lowered. In the Fig. ^, we can 
see that the curves for T = 104K and T = 52A" cross 
around kpr — 3.2. The temperature where K(r, T) has 
the maximum value is a function of separation r and has 



lower value with larger separation r. The Knight shift 
has different r dependence whether temperature is above 
the Kondo temperature or below the Kondo tempep.- 
ture as mentioned by E.S. S0rensen and I. Afheckc^rEa. 
The Knight shift converges faster at higher temperatures 
(above the Kondo temperature) and it has longer range 
at lower temperatures. The Knight shifts of Yb ions at 
fixed angles show similar behavior. 

For a fixed Kondo temperature To and distance r, 
the knight shift K(T) shows a linear relation with the 
magnetic susceptibility x(^) a ^ high temperature and it 
starts to deviate from the linear relation and shows an 



anomaly when temperature is lower than T, 
x(T) reaches its maximum value. Fig. 
Knight shift at 9 = 0, T = 380K as 



where 
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shows the 
a function of 

the magnetic susceptibility for different separations, with 
temperature as an implicit variable for Ce ions. These 
curves show a linear K vs. x relation at high temper- 
atures (in this figure, high temperatures correspond to 
the lower left corner) and show non-linear K — vs. — x 
relations at low temperatures. The anomaly, i.e. the 
magnitude of the non-linearity, diminishes as the sep- 
aration r is decreased. In Fig. [l^, the Knight shift is 
investigated at kpr = 3.3, and 9 — for different Kondo 
temperatures. As in Fig. O, temperature is an implicit 
variable and high temperatures correspond to the lower 
left corner. The magnitude of deviation from the linear 
relation is decreased when the Kondo temperature is re- 
duced. For CeSn3, kp^k ~ 25 where £r- = hvp/ksTo 
is the Kondo screening length with the Fermi velocity 
vp. So our calculation is done well inside the conduction 
electron screening spin cloud. Very similar results are 
obtained for Yb compounds. n 
This results qualitatively confirm the Ishii's arguments 
that when the radius r is bigger than Kondo screening 
length = hvp/kgTo, where vp is the Fermi veloc- 
ity, the anomalous conduction spin density cloud sets in. 
Far inside this radius, conventional temperature indepen- 
dent RKKY oscillations dominate of the kind observed 
by Boyce and Slichter. Outside the screening length, 
at T = 0, the anomalous term will dominate also with 
an RKKY form but an amplitude of order D, the con- 
duction bandwidth, compared with D(N(0)J) 2 for the 
Ruderman-Kittel term, where N(0) is the conduction 
electron density of states at the Fermi energy and J the 
conduction electron- local moment exchange coupling. 
Ishii did not calculate the explicit temperature depen- 
dence of this structure, but did anticipate that it would 
vanish above the Kondo scale. Scaling analysis confirmed 
the asymptotic factorization oLthe Knight shift for short 
distance and low temperatures. A possible understand- 
ing of the Boyce and Slichter results, then, is that the 
Cu nuclei they sampled were at distances r << £k from 
the Fe ions. But, in our calculation, the anomaly is still 
present and the magnitude is surprisingly large for short 
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distance provided To is large. A heuristic basis for under- 
standing this is the two-pole approximation, as discussed 
in subsection III E . 



B. CeSna 

The compound CeSn3 has the f.c.c. AuCii3 crystal 
structure. The local symmetry at every cerium site is 
cubic. With decreasing temperature, the magnetic sus- 
ceptibility x(T) of CeSn 3 shows first typical Curie- Weiss- 
like behavior, followed by a maximum aL.T max « 40K 
and tending to a constant value at T = Q3l1 CeSn3 has 
positive amplitude of Knight shift for the 119 Sn nuclei. 
However, our calculated amplitude of the Knight shift 
before scaling it to x is actually negative. This implies 
that the fit is sensible only if the assumed contact cou- 
pling between conduction and nuclear spins is negative. 
This actually makes sense because the Sn nucleus should 
dominantly couple through core polarization, which pro- 
duces a negative effective contact coupling. The NMR 
Knight shift K(T) and x(T) are related linearly above 
T m ax- But the Knight shift doesn't follow the magnetic 
susceptibility below T max and has weaker temperature 
dependence . The 119 Sn NMR Knight shift has a positive 
anomaly, i.e. it has larger magnitude than that of mag- 
netic susceptibility below the T max . The positive muon 
spin rotation(/isr) measurement shows that the positive 
muon Knight shift also exhibits an anomaly below T max , 
but it has a different sign with respect to the susceptibil- 
ity as compared to the NMR Knight shiftQ. 



of the thermodynamics in many cases. 

Because Sn is on the face of each unit cell, we also aver- 
age the Knight shift over the reference frame determined 
by whether 119 Sn is on the xy plane, or the xz plane(for 
H || z). Note that yz is equivalent to xz. Like the 
systematic calculation, a Van-Vleck term is included for 
both magnetic susceptibility x{T) an d the Knight shift 
K(T) calculation. The Knight shift is scaled by an inter- 
mediate range temperature match to the susceptibility, 
because the experimental data is measured only up to 
room temperatures. 

As shown in Fig. |l3| the results agree well with the ex- 
periments in spite of the oversimplified conduction elec- 
tron band structure. Experimentally the 119 Sn NMR 
Knight shift has a positive sign and is linearly related 
to Ce magnetic susceptibility x{T) at high temperatures. 
K{T) shows an anomalous deviation from Ce x(T) at low 
temperatures and the magnitude is bigger than x- At 
high temperatures above 4QQK : K(T) is reduced slowly 
with increasing temperature and shows deviation from X- 
Experiment measures the Knight shift only up to 300JT, 
so we are not sure whether this behavior appears in ex- 
periment. The magnitude of the anomalous contribution 
goes down with distance from the Sn nucleus and the the- 
oretical data at a fixed distance which most closely match 
those of experiment are taken fcj?r = 2.1 and 6 = 0. Note 
that this distance is an order of magnitude smaller than 
the Kondo screening length £r- = Tivf I^-bTq ~ 2hjkp- 



D. Lattice Coherence Effects 



C. NMR Knight Shift 

To study the Knight shift for CeSn3, first the hy- 
bridization width r for CeSn3 is decided by comparing 
the calculated magnetic susceptibility and experimental 
magnetic susceptibility, assuming a conduction electron 
density of states half width D = 3eV. As in the sys- 
tematic calculation, a fixed value for parameter e/5/2 = 
-2eV = -0.6667D is assumed, and Aso = 0.29eV = 
0.096667-D is used. Then D is varied a little to fine tune 
the results. The Knight shift is calculated with the fixed 
Kondo temperature (specified by F = 0.050666D) and 
conduction band width D — 2.655eV. The best fit Kondo 
temperature which is given by Eq. ( |67|) is To = 0.0123, 
in units of D i.e. T = 380K. 

For the Knight shift, an incoherent lattice sum is car- 
ried over several hundred shells of surrounding Ce ions 
about a given 119 Sn nucleus. We assume the Knight shift 
contribution of each ion to be described by this single im- 
purity model, known to be a good approximation at high 
temperatures where the ions are incoherent with one an- 
other, and known to provide a very accurate description 



We study lattice coherence effects within a local 
approximation^ = 00 expansion) to the lattice model. 
The conduction electron propagator in Fig. |5| becomes 
dressed and the conduction electron self energy is calcu- 
lated using the average T-matrix approximation, which 
is exact for a Lorentzian density of states within the lo- 
cal approximatioroEJ, and otherwise corresponds to the 
first iteration of the d = 00 self consistency. The same 
parameters for the incoherent Knight shift calculation are 
used for the coherent lattice calculation. 

The calculated results fit experimental data well in the 
region of temperatures between 10QK and 25QK where 
the experimental data exist. At low temperature, the 
Knight shift with coherent lattice effects also shows an 
anomaly, with a little discrepancy with the experiment. 
At higher temperatures it shows tails which has bigger 
values than the experimental susceptibility (see Fig. fl7| ). 
If we fit the results to this high temperature magnetic 
susceptibility, then the Knight shift has bigger values at 
the maximum temperature and shows sharp changes with 
temperature changes. We are not sure which is the best 
way to fit the result with the experiment. The calculated 
magnitudes (before the scaling to fit the experimental re- 
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suits) of both incoherent and coherent Knight shift have 
similar values at high temperatures where coherent lat- 
tice effects are small. Explicit values are shown in Ta- 



ble in. 



The average magnitude and amplitude of the oscilla- 
tion of the Knight shift is decreased from the that of 
the incoherent lattice sum; we believe because of the 
damping effects brought in by the imaginary part of 
conduction electron self energy. To test this idea, we 
added a phcnomcnological constant damping to the in- 
coherent lattice sum. Fig. [l4| shows the result for small 
damping , — ImS/ = 0.01D and Fig. [l5| shows the re- 
sults for large damping, —ImS/ = ID. In these figures 
K ( T ) = Ek F r<k Fro K(r, T). This study shows that the 
amplitude of oscillation is reduced when the damping is 
increased. And the amplitude of the Knight shift con- 
verges the faster when the damping is the bigger when 
distance is increased. An impurity at large distance 
doesn't contribute to the total Knight shift because of 
the large damping, i.e. short life time of the conduction 
electrons. Also Fig. [l6] shows the Knight shift calculated 
with the coherent lattice sum. It shows a behavior inter- 
mediate between small and large damping calculations. 



E. ^sr Knight Shift 

Positive muons, stopped in a solid, come to rest at 
interstitial sites where the muon spin performs a Lar- 
mor precession in the local magnetic field. The muon 
Knight shift is then a measure of the local magnetic sus- 
ceptibility. For CeSn3, from volume considerations it is 
most likely that the muon preferentially occupies the oc- 
tahedral interstices of the AuCu3 structure (as it does in 
metals with the closely related fee structure), rather than 
the tetrahedral sites. There are, however, two inequiva- 
lent octahedral sites, one at the center of the cubic unit 
cell and the other at the middle of the Ce atoms which 
has a non-cubic symmetry. So in principle two resonances 
are expected. In the experiment by Wehr and et alt!3, 
only one resonance was observed because either the muon 
performs a site average by fast diffusion or the frequency 
difference is small with respect to the apparent width of 
the signal which is given by the muon lifetime and the 
intrinsic width 7. 

In the Curie- Weiss regime of CeSn3 above 200K, the 
temperature dependence of the positive muon Knight 
shift is linearly related to the bulk magnetic susceptibil- 
ity. In the intermediate valence regime of CeSn3 the lo- 
cal magnetization as experienced by the muon decreases 
more strongly below 200K than the magnetization of the 
4/ state as deduced from the bulk susceptibility. This 
behavior reflects either a modification of the transferred 
hyperfine fields between the 4/ moments and the muon 
or signals the influence of an additional negative Knight 



shift contribution which was absent or small in the high 
temperature ranges. Comparing with the NMR Knight 
shift, the /zsr Knight shift has a maximum at higher 
temperature and the sign of anomaly is opposite. This 
anomalous reduction of positive muon Knight shift might 
be regarded as an indication for an_additional negative d- 
electron Knight shift contribution!^. In terms of a band 
picture, the increase of d character at the Fermi level can 
be understood as a 4/ — 5d hybridization effect. It is sup- 
ported by de Haas-van Alphen measurements!!! Positive 
muons sitting between Ce atoms should be particularly 
sensitive to variations in the 5c? states from the symmetry 
of these orbitals. 

We have calculated the fj, + Knight shift both possible 
muon sites. Both incoherent and coherent lattice sums 
are carried out over several hundred shells of Ce atoms. 
All other parameters for the calculations are same as for 
the 119 Sn NMR Knight shift calculations. Fig. 18 shows 
the results assuming that muon sits at the center of the 
cubic unit cell. At high temperatures, the calculated [isr 
Knight shift agrees well with the experimental data and 
shows a linear relation with the bulk magnetic suscep- 
tibility. Both incoherent and coherent lattice sum stud- 
ies give the correct magnitude for the low temperature 
Knight shift anomaly but the wrong sign. Also the max- 
ima occur at lower temperature than the data and the 
magnitude is bigger than the experiment. 

Results are shown in Fig. [l9] for the case where that 
muon sits at the middle of the Ce-Ce band axis. In this 
case, results with the incoherent lattice sum show sim- 
ilar behavior to the previous calculations (assuming the 
muon sits at the center of the unit cell). But the result 
of the calculation shows very interesting behavior even 
though it does not agree with experiment. The Knight 
shift starts to deviate from the linear relation at higher 
temperature and the magnitude of the anomaly is larger 
than the experimental result. Note that the sign of the 
anomaly is agrees with experiment. There is a possibility 
to fit the experimental data by averaging the fisi Knight 
shift from both positions. Because there is no experimen- 
tal data which gives information about the fractional site 
occupancy we just averaged two Knight shifts with sev- 
eral fractional occupancy ratio / which is defined as 



^total = /-^center + i 1 ~ f) ^bond- 



axis ' 



(68) 



Fig. |2Cj shows the result for / = 2/3. Our calculation 
misses the maximum point, but the low temperature 
anomaly agrees with the experimental result. 

Also, if the muon is not situated at a site of cubic 
symmetry, dipolar fields from the induced local moments 
may give the dominant contribution to the positive muon 
Knight shift. The magnetic dipole which is inversely 
related to the mass, can give comparable contribution 
to the positive mupn Knight shift contrast to the 119 Sn 
NMR Knight shiftEa. The direct dipolar interaction en- 
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ergy of two magnetic dipoles fhi and fh-z, separated by r 
is given 

U = -g [mi ■ m 2 — 3(mi • f)(m 2 • f)] (69) 

The dipolar part of the Knight shift tracks directly the 
atomic susceptibility of the nearest neighbor f-electron 
ions while the contact hyperfine contribution is what we 
are interested in. In our calculation the dipolar field effect 
is of course exactly cancelled out when we average over 
the reference frame. 



incoherent lattice sum is carried out for over 8000 shells 
of atoms, larger than for CeSn 3 . Because of the com- 
plicated crystal structure, each shell only includes a few 
atoms. 

Calculated results are shown in Fig. |2l]. Both the mag- 
netic susceptibility and the Knight shift agree well with 
experimental results and can explain both the magnitude 
and sign of Knight shift anomaly in spite of the oversim- 
plified conduction electron band structure. We note that 
the sign of the anomaly is opposite to that of Ce in this 
case, and indeed we find that these contributions go in 
opposite directions numerically. 



F. YbCuAl 

The ternary inter-metallic compound-YbCuAl has the 
hexagonal Fe2P type crystal structureOEj, in which each 
Yb atom has the same local environment. 

At low temperatures, the magnetic susceptibility x(T) 
has a large constant value(x(0) = 25.5 x 10 _3 e. m. u. I mole 
Yb atoms) and a maximum value at T max w i^S^-i There 
is a Curie-Weiss like behavior above T maK llIrLll. 27 A1 
NMR shift data were obtained from derivative spectra 
of the central (1/2 <-> -1/2) NMR transition. Above 
X and K track each other, as expected if only one 
mechanism is appreciably temperature dependent. Here 
the Yb magnetization is the obvious candidate for the 
temperature-dependent contributions to both x(T) and 
K(T). This linear K — vs. — \ relation has been used to 
determine the relative scales of the x(^) and shift coor- 
dinate axes in Fig. ^l|. The linear relation breaks down 
below T max . 

The 27 A1 NMR Knight shift has negative sign and the 
absolute magnitude of the low temperature Knight shift 
is smaller than the magnetic susceptibility, opposite to 
the case of CeSn3. For YbCuAl, the ground state energy 
for f state is taken as e/7/2 = —leV = —0.3331? and 
A so — 1.3eV — 0.4333D. Because of the large value 
of A so , we can neglect the interaction between j = 7/2 
ground state and j = 5/2 excited state. Without this 
Van-Vleck term, we can estimate the conduction electron 
band half width D using the zero temperature magnetic 
susceptibility x(0) value. 



x(o) = 



gjjij + l)u 2 B N 
3k B T 



(70) 



With x(0) = 25.5x 10~ 3 , we get k B T Q = 8.7974x 10" 3 eF. 
Also with the relation 



(71) 



/ //< 



where D 



the 



3eV, and (T /D) the = 0.00184, we get 
D eX p = 4.7812eV. But, we get the best magnetic sus- 
ceptibility fit with D = 3.5eV and T = 0.02033D. The 



Cl( 

al. 



G. Y0.sU0.2Pd3 

The Yi-^U^Pda 3 system has aroused great interest, 
following the discovery of non- Fermi licmicLhehaviour for 
uranium concetrations around x = 0.2c3'E 2 iL 3 J. This de- 
screpancy has been interpreted as possihly-|arising from a 
two channel quadrupolar Kondo efTectE3'EZl or from criti- 
cal effects of a new kind of second-order phase transition 
at zero temperatureEl We will mainly discuss the com- 
position x = 0.2. 

Y0.8U0.2Pd3 has a cubic AuCu3 crystal structure. The 
j = 4 ground state of U 4+ is split to Ti singlet, T 3 
nonmagnetic doublet, and T4 and T$ magnetic triplet. 
Knowledge of the crystal field ground state is a crucial 
test for the validity of the quadrupolar Kondo model and 
there are several neutron scattering qYflPr.iments to de- 



thc ground state of Yo. 8 Uo.2Pd 3 BEl& Mook et. 
!3 interpreted their results in terms of a crystal field 
level scheme with a T3 doublet ground state and T§ and 
T4 excited triplet states at 5 and IdmeV, respectively 
and thus support the two channel quadrupolar Kondo 
effect interpertation. In this case, x(^) originates in the 
Van-Vleck susceptibility associated with transition from 
a T3 nonmagnetic ground doublet into excited-state Ts 
and T4 triplets The T3 is described by a quadrupolar 
pseudospin. This couples to pseudospin variables of a 
conduction Tg quartet in time-reversed channels with the 
antifcrromagnetic pseudospin coupling mediated by vir- 
tual charge fluctuations to magnetic doublets in excited- 
state configuratiops. 

McEwen et alE2 saw a peak of magnetic origin at 
36meV and another peak around AmeV and explained 
this with the transition between T3 ground state with 
excited states T4, Ti and Y§. They couldn't find a peak 
at 16meV. _ . 

Dai et alSfH reported a T5 magnetic ground state with 
polarized inelastic neutron scattering experiment and T3 
and T4 excited states at 5, 39 meV. We note, however, in 
contrast to this interpretation, that there is no quasielas- 
tic scattering around Tioj — 3.6meV(~ A2K = /c;,Xb), 
expected for a conventional magnetic Kondo effect. 
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For x = 0.1 and 0.2, a breakdown of the expected 
linearity between the NMR Knight shift-and the bulk 
susceptibility \ is found below ~ 50^19. The mag- 
netic susceptibility exhibited an upturn at low temper- 
atures(Curie tail), indicating the presence of magnetic 
impurities(see Fig.g2h_. The impurity magnetic suscepti- 
bility_was subtractedED. The temperature dependence of 
x(T)E3 suggests that the mechanism for the two channel 
behavior is the quadrupolar Kondo effec1c3. 



x(o) V 



(72) 



In our study, we consider f 2 and / 3 configurations for 
U ions and only the Hund's rule ground states, i.e. / 2 , 
j = 4 and / 3 , j = 9/2 spin-orbit states are kept for the 
calculation, j = 4 states is split to ground doublet, 
T5, T4 and Ti excited states and j = 9/2 multiplet is 
split to Tg doublet and two Tg quartets. The conduction 
electron band width D is assumed to be 3eV. All param- 
eter values are listed in Table |ll] in the unit of D. The 
incoherent lattice sum is carried out over three hundred 
shells with impurity configuration averaging. 

Fig. ^2] shows both the experimental and calculated 
magnetic susceptibility and 89 Y NMR Knight shift. In 
our calculation, both the bulk magnetic susceptibility 
and the Knight shift become constant when the tempera- 
ture goes to zero and thus a Knight shift anomaly doesn't 
arise. Our calculated magnetic susceptibility saturates 
when the temperature goes to zero and doesn't show 
the low temperature singularity like experiment. This 
may arise from the numerical calculation or from intrin- 
sic properties such that real ground state may be mag- 
netic as discussed earlier. A separate possibility is that 
the weak admixture of excited / 3 magnetic states con- 
tributes a weakly log divergent contribution to x; this 
possibility may be explored elsewhere. Regardless, the 
calculated Knight shift agrees well with the experimen- 
tal data in magnitude and temperature dependence. 



V. SUMMARY 

In this paper, we have calculated the magnetic sus- 
ceptibility and the Knight shift for the heavy electron 
materials within the infinite-?/ single impurity Anderson 
model using NCA method. 

In our calculations we can explain that the Knight shift 
anomaly in heavy electron materials with the simplified 
single impurity Kondo effect. There exists a Knight shift 
anomaly at short distance r < £k, with amplitude pro- 
portional to Tq/T. 

Our calculations show generally good agreement with 
experimental results in spite of the oversimplified band 
structure. Especially, the short distance Knight shift de- 
pends on the detailed structure on conduction electron 



band and our calculation shows large contributions from 
the short distance Knight shift. For future work, we can 
include more realistic conduction electron band structure 
which can be calculated with LMTO (linearized Muffin- 
Tin Orbital) method. 



ACKNOWLEDGMENTS 

This research was supported by a grant from the U.S. 
Dept. of Energy, Office of Basic Energy Science Division 
of Material Research. We thank H.R. Krishna-murthy, 
D.E. MacLaughlin, M. Steiner, and J.W. Wilkins for use- 
ful conversations on this work and related topics. 



APPENDIX A: U 4+ IONS 

For the U compound , the crystal field effect (CEF), 
which lift the angular momentum degeneracy of U ions 
and their spin-orbit multiplet decomposes into irre- 
ducible representation of the cubic field need to be in- 
cluded. The distinction from the CeSn3 and YbCuAl 
cases is that the apparent crystal field splitting Acef ~ 
T for CeSn 3 and YbCuAl. The f 2 ,j = 4 Hund's rule 
ground state of the U ion is split to a T3 nonmag- 
netic doublet, r§ and T4 magnetic triplet and Ti sin- 
glet statesc 3 ]. The f 3 ,j = 9/2 spin-orbit multiplet is 
split to Tj doublet and two Tg quartets. For an explicit 
derivation, see Appendix [BLand the articles by K.R. L 
M.J. Leask and W.P. Wolfel and by M.T. Hutching 
The eigenstates of 1^ CEF states for j ^-A^znd j = 9/2 

]'L2L In cubic 



multiplets are in Table IV and Table 



symmetry, the coefficients of CEF states depend upon 
the parameter x which is fixed by the ratio of the fourth 
and sixth degree terms in a short distance expansion of 
the cubic field in the Hamiltonian of the crystal electric 
field, and upon the parameter W which is an overall scale 
factor fixed by the crystal field strength. In this calcula- 
tion, we use x 3 = 0.3693 and W 3 = 2.746 x 10~ 4 to have 
j = 9/2 Tq for the ground state of the f 3 configuration 
and x 2 = -0.648 and W 2 = -3.95 x 10~ 4 to have j = 4 
T3 for the ground state for f 2 configuration. For further 
details, see Appendix [b|. The choice of the overall phase 
is arbitrary in defining the CEF eigenstates. 

A brief descriptions of the different irrep labels of the 
cubic group is as follows: 

(1) Ti and T2 are orbital singlets. 

(2) F3 is an orbital (non-magnetic, or non-Kramers') dou- 
blet and usually labeled by ±. 

(3) T4 and are magnetic triplets and labeled by 0, ±1. 

(4) Tg and T7 are magnetic Kramers' doublets and la- 
beled with pseudo-spin f or J,. That is, the Tg and 
CEF states are similar to the j = 1/2 angular momen- 
tum manifold. 
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(5) T 8 is a magnetic quartet (T 8 = T 3 ® Ty). and labeled 

by ± T / 4- 

In the Anderson model picture, the conduction elec- 
trons can hop on and off the atomic orbitals at the impu- 
rity site, the / = 3 conduction electron partial waves are 
most strongly coupled to the / electrons in U ions (for the 
isotropic hybridization, only the I — 3 components can 
hybridize with the / orbitals.). In the presence of the 
spin-orbit coupling, the I — 3 conduction electron state 
splits into the j = 5/2, 7/2 conduction electron states. 
These j multiplets further split into the CEF irreducible 
representations, Tq, T7 doublets and Tg quartet in crys- 
tal environment. The CEF eigenstates for j = 5/2, and 
j = 7/2 multiplets are listed in Table and Table [yTi . 

The hybridization Hamiltonian of our model U com- 
pound in the absence of the CEF is given by 



APPENDIX B: CRYSTAL ELECTRIC FIELD 
EFFECT 

In the crystal lattice, magnetic ions feel an electrostatic 
field produced by the neighboring ions. This crystal field 
lifts the degeneracy of the angular momentum of the mag- 
netic ions. The most common method to calculate the 
effect of thexrystal electric field is the operator equivalent 
tcchniquesE3 which exploits the Wigner-Eckart theorem 
to replace the electrostatic potential terms in the Hamil- 
tonian by operators in the angular momentum space of 
the ground multiplet. It depends on the symmetry of 
the crystal and the orbital angular momentum j of the 
magnetic electrons— .The most general Hamiltonian with 
cubic symmetryS'O can be written as 



H 



cef 



B 4 (O 4 



50l) 



B 6 (O 6 



21Q*) 



(Bl) 



H 



of 



E 

kj c m c 



E 



Hj c m c ;f 2 j 2 = 4m 2 |/ 3 j 3 = 5^3^, 



|/ 2 i2 = 4m 2 )(/ 3 i 3 = -m 3 \+h.c.]. 

Here the reduced matrix elements are 

9 

A(j c m c ;/ 2 j 2 = 4m 2 |/ 3 j 3 = -m 3 ) 

9 

= (fj2 = 4TO 2 |/ icm J/ 3 j 3 = 2™3> 

= K(j c j 2 1 j 3 ) (j c m c ; j 2 m 2 \j 3 m 3 ) 



(Al) 
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ol 



354 - [30 J(J 
-6J(J+ 1) + 3J 2 (J 



1)-25]J2 

2 / 7 1 1 \2 



If 



oi 

ol 



V 



+ J 4 -) 



231J Z 6 - 105[3J(J + 1) - 7] J\ 



(A2) 



+ [105J 2 (J- 
-5J 3 (J + 1) 3 

0\ = \\\\Jl- J(J 



l) 2 - 525J(J+ r 



5J 3 (J + 
+ 294]J 2 



40J 2 (J + l) 2 - 60 J(J+ 1) 
l)-38](4 



where j c = 5/2,7/2. The Wigner-Eckart theorem is 
used for the last line and the prefactor K(j c j 2 \j 3 ) is 
the fractional parentage coefficient. In our calculation 
where j 2 = 4 and j 3 = 9/2, K(j c j 2 \h) = \/V§ for 
both j c = 5/2 and j c = 7/2. {j c m c \ j 2 m 2 \j 3 m 3 ) are the 
Clebsch-Gordan coefficients. If there is a crystal electric 
field, the Anderson hybridization between the CEF states 
Ts needs to be re-evaluated in the basis appropriate to 
the crystal field split states. We have 



Jt) 

+ -U% + Jl)[lU 2 - J(J+ 1) - 38] 



(B2) 



! ''< r= E 

fej c r c a c Tar/3 



Hcf = 



E t V"C 



kj c V C 



A( Jc r c a c ;4ra|-r/3) 



The coefficients B4 and B e are the factors which de- 
termine the scale of the crystal field splittings. In a 
simple point charge model, they are linear functions of 
< r 4 > and < r 6 >, the mean fourth and sixth powers of 
the radii of the magnetic electrons, and thus depend on 
the detailed nature of the magnetic ion wave functions. 
We treat these as phenomenological parameters because 
these are very difficult to calculate quantitatively. 

Following ReE3, we rewrite the Hamiltonian as 



\fHTa)(f-T[3\ 



h.c] 



Hcef = B A F(A) 



Q 4 
F(4) 



B 6 F(6) 



Q 6 

F(6)' 



(B3) 



E 

fcj c r c Q 



Ei 

Tar/3 



V k K(j c 



4Ta|-r/J> 



where 

4 = {(oi 



50|], and 6 = [0\ 



2101} 



and F(A) and 



-kj c r c a, 



\fHra)(f-r/3\ + h.c.}. 



(A3) 



The reduced matrix A(_j c r c a c ; 4ra||r/3) is implicitly de- 
fined in above equation. All the possible selection rules 



F(6) is the common factors to all the matrix elements of 
fourth and sixth degree terms. In order to cover all the 
possible values of the ratio between the fourth and sixth 
degree terms, we put 



for the hybridization are listed in Table VIH 



B 4 F(A) 
B e F(6) 



Wx 
W(l 



M) 



(B4) 
(B5) 



where — 1 < x < 
It follows that 
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B 6 



F(6) 



(B6) 



so that B^/Bq 
x = ±1. 



for x — and B^/Bq = ±oo for 



Rewriting Eq. (B3) we have 
4 



H 



cef 



W[x{- 



(1 



W)(-^t)1- 



(B7) 



Here, the eigenvalue is related the crystal electric energy 
by the scale factor W . 

In our calculations, we need the relations between pa- 
rameters X2 and W2 for f 2 configuration and X3 and W3 
for f 3 configuration of U ions. To get j = 4 T3 and 
j = 9/2 T 6 doublets as the ground states of / 2 and f 3 
configuration, we choose x-i < and X3 > and W% < 0. 
W2 is decided by the energy splitting between T 3 and T 5 
states, A35. 



T^ 2 = -A 35 /(84+ 114x 2 ) 



(B8) 



,/From the Eq. flB4|), Eq. ( p35[) and Eq_(|B6|), and factors 
which is given in Table 1 of reference l£j, we get the rela- 
tions of parameters between f 2 and f 3 configurations. 



3.143i 3 (l - M) = s 2 (l - \x 3 \) 
W 3 = 0.39631^2 — 

X3 



(B9) 
(BIO) 



^From the above equations we get x 3 = 0.3693 and W3 ■ 
2.746 x 10~ 4 for x 2 = -0.648 and W 2 = -3.95 x 10~ 4 



APPENDIX C: EVALUATION OF KNIGHT 
SHIFT K SUM 

In this appendix, the form of -T(e&, ey) expressed in Eq. 
( |36| ) is derived. We define 

'.//•"/..'/.<i ' /• ! l/.n/ :i iA' , rj-i 



h(kr)j 3 (k'r)^ £ G (fc CT , iw)G (*V,tw) 



Let's first do the v summation. 
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(C3) 



where /(p)[iV(C)] = (e /3p + l)" 1 [(e /3C - l)" 1 ]. In this Eq. 
(|C3]), only {p — ( — iuj) term is depend on lu. Let's do the 
u sum. We get 
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Now taking the A — ► —00 projection 
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Also with the relations 
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The above equation, can be rewritten 
ja{kr)j3[k'r) 



(C7) 



where uj\v\ = [2n + l)m / f3[2rnri / 0\ and the usual 
fermion(boson) Matsubara frequencies. Go is the bare 
conduction electron propagator, Gj m is a pseudofermion 
propagator for spin-orbit multiplet j, and D is a pseudo- 
boson propagator. For the coherent calculation the bare 
conduction propagator is replaced by the dressed electron 
propagator. This coherence effect will be considered in 
Appendix |]. 



x (f(e k )h djl (e k ) + (1 - f{e k ))h r y (e fc )) (C8) 



where 
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Then Eq. (CI) can be rewritten as 
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J it iv — C J ir i 
dp' A jlm , (pQ 
7r + — p' 



^jm(p) 



((J + v ) — p 



J_ /' d P /1 (- 
24/. 



- / ^Aj m ) (p)i? e G /m ,(p)i? e ^(p- £fe ). 

■7 7T J 



(C2) 



(CIO) 
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APPENDIX D: ANALYTIC CALCULATIONS OF 
INNER K 1 INTEGRAL IN THE KNIGHT SHIFT 

We can calculate the inner kl integral analytically. 
This is given by 



k 2 - k' 2 



(Dl) 



where j n is the spherical Bessel function. The above 
integration is broken into 4 terms which may be expressed 
in terms of Sine and Cosine integrals. 
The first term is 



dk' 



k' 2 cos k'r 
k 2 — k' 2 k'r 



1 7T 

= — \Ci(kr) cos kr + si(kr) sin kr H sin kr] 

r 2 

where the Sine and Cosine integrals are 
sini 



Si(z) 



t 



-dt 



si(z) = Si(z) 



Ci(z) = 7 + lnz + 



cos t — \ 



t 



-dt. 



(D2) 

(D3) 
(D4) 
(D5) 



For more explicit-formalism of Sine and Cosine integrals 
see the Referenced. The second term in our integration 



is 



dk'-. 



-6 sin k'r 



• - - k' 2 (k'r) 2 
■ [Ci(kr) sin kr — si(kr) cos kr 

The third term is 

k 12 — 15 cos A/ r 



1.2 

o K 
-6 

kr 2 



■cos fer]. (D6) 



dk' 

k 2 -k' 2 (k'r) 3 

-15 



k 2 r 3 
ir 



Ci(0) + Ci(kr) cos kr + si(kr) sin kr 
sin At]. (D7) 



Ci(0) is divergent logarithmically, but is exactly canceled 

by a term from next integral. 

Our fourth term in the integral is given by 



dk'- 



k' 2 15 sin A:' 



„ k 2 -k' 2 (k'r) 4 

15 7T 

■: \Ci(kr) sin kr — si(kr) cos kr H cos kr] 

k 6 r 2 



15 

k 2 r 3 



[l-tft(O)]. 



(D8) 



Note the explicit cancellation of Ci(0) between Eq. (D7) 
and Eq. (Q. Putting Eq.(p|)-Eq.® together, we 



have our inner momentum integral in our Knight shift 
expression define in Eq. (Dl) given by 



k' z dk' 



k 2 - k' 2 



1 7T 

= — [Ci(kr) cos kr + si(kr) sin kr H — sin kr] 
r 2 

6 . 7T 

\Ci(kr) sin fcr — si(kr) cos /cr H cos kr] 

kr 2 2 

15 7T 

5— r [Ci (kr) cos kr + si(kr) sin fcr H sin fcr — 1] 

k 2 r A 2 

15 7T 

+ [Ci(kr) sinfcr — si(kr) cos fcr + — cos fcr]. (D9) 
fc d r 4 2 



APPENDIX E: ANGULAR DEPENDENCE OF 
THE KNIGHT SHIFT 

In this appendix the derivation of fjj> (a, 8), the angu- 
lar dependence of the Knight shift, will be discussed. 

First, we need to calculate the expectation value of the 
magnetic moment operator in the z direction 

(jrrij\J z + S z \j'm. r ) = rrijSjf + (jmj\S z \j'm. r ). (El) 

Then for j — j' , by Wigner-Eckart theorem gives 

(jrrij]J z + S z \jrrij>) = 'i'„. „. ..'/,'";• (E2) 

where gj is the Lande g-factor and g 5 / 2 = 6/7 and 57/2 = 
8/7 for I — 3 and s = 1/2 which is the case for 4/ 1 and 
4/ 13 configuration. And for j ^ j' 

(jmj\J z +S z ]j'm j >) = (jm^S^j'mf) 

= S mj m jl a(jmj\lmy, sa){lmy, sa\j'mj>) (E3) 

mi a 

where a = ±1/2. For I = 3, s = 1/2, j = 5/2 and 
f = 7/2 

(■^mj\S z \-m,j>) = S mjm .,^2 a(-m J \3m 3 ;sa)(3m 3 ;sa\ 



"TO, 



7713a 



49 - 4to^ 



14 "mjOTji 



(E4) 



given 



(jim - a; ^a|i, m) = ( ^^"'f ) 1/2 for j = j x + 1/2 (E5) 
2 2ji + 1 

-^ JlTm/2 ) 1/2 forj=^-l/2 (E6) 



2ji + 1 



where 



(3to 3 ; \ ± i|Im> = (^l^) 1 / 2 for j = 7/2 (E7) 
(3to 3 ; \ ± = T(^|^) 1/2 for j = 5/2. (E8) 
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Then 



{\ m j\Sz\\ m ]> 



7 



for rrij = ±1/2 
for rrij = ±3/2 



= -^— iovmj =±5/2 (E9) 

In the Knight shift calculation the angular dependence 
comes from the Zeeman term in the Hamiltonian. 

H z = -H B H Z {L Z + 2S Z ) = -fM B H z (J z + S z ), (E10) 

where external magnetic field H = H z z. 

Now let's calculate the angular part of the Knight shift. 
To do the lattice sum, we have to consider the difference 
of the field axis and the bond direction r which con- 
nects the nucleus or muon to a given / ion. Let the 
angle between the field axis and r is a, and the an- 
gle between the z axis and bond axis is 9. First when 
the field is along the z direction, i.e. 9 = a and nu- 
clear spin I = 1/2, the angular momentum operator 
J z , which is quantized in the bond direction becomes 
cos9J z — sin9J x — cos9J z — sin9(J + + J-)/2 in the 
new reference frame when the material has the cubic 
symmetry. Also, the nuclear spin operator a z becomes 
cos 9a z — sin 9a x = cos 9a z — sin 9(a + + er_)/2. Then the 
surviving terms in the Knight shift calculation are 

cos 2 9a z (J z + S z ) + sin 2 9a x (J x + S x ) = 
sin 2 9 

cos 2 9a z {J z + S z ) + -^-(a+(J_ + S-) + <7_(J+5+)). 

(Ell) 

Then the total angular part is 
f jA 0) = cos 2 9f h , + +fc + f$ + ®). 



(E12) 



Where /?■, is 

fjj'( 9 )= E ti m j\ J z + s z\? m 3')<a(? m j'\r)(r\jri 

rrij rrij / a 

= E a U m j\ J z + S z \j'm r ) (j'm r |3m 3 ; ]-a) 



rrij ams 



x(3m 3 ; ■za\jm j )\Y 3ma (r)\' 



(E13) 



where a = ±1/2 and Yi m (f) is the spherical harmonics. 
Then for l)j = f 

fjj( e ) = E a(jm J \J z + S* 2 |jTO^-)(jm^|3m 3 ; ia) 

rrij rrij arris 



x(3m 3 ;-a\jmj)\Y 3m3 (f)\' 



(E14) 



With the above equation and Eq. (E7) and Eq. (E8) 

/l = W = ^-(l-4sin 2 0) (E15) 

2 2 ZiiTT 

/lz(0) = ^(l + 3sin 2 #). (E16) 

2 2 Y7T 

2)j i.e. j = 5/2 and/ = 7/2 

This term gives the van Vleck Knight shift 

$\l^ = E a (j m j\ J z + S z \3'mj'){j'rnj<\3m3-,-a) 

rrij arris 



x(3m 3 ; ^a\jmj)\Y 3m3 (f)\' 



Utt 



(2 -sin 2 9). 



(E17) 



Now let's derive f^p, the angular part of Knight shift 
from the second term of Eq. ( EllD , where j = j' 



rrij rrij/ a/3 

E 0'' m jl J ±li'' m J')0''™i'|3m 3 ; -a) 

rrij rriji msm'^ a/3 

x(3m 3 ; ^a\a^\3m 3 ; ^/3)(3m 3 ; -(3\jmj) 
xY 3 * m3 (?)Y 3m , 3 (r) (E18) 

with 

(j,m j \J±\j',m jr ) = ± mj)(j T rrij + l)8 jr 5 mjim] ^i 
and 

(3m 3 ; ^a\a^ t p\3m 3 ; = \Jj~ P((3 T l)6 a p Tl Sm 3 m' 3 

The above equation gives 
fifW = E ± m M j T ^Tl)0'^Tl|3m^; ~ ±~) 



1 1 



x(3m 3 :-±-|im i )|Y 3m . T i(f)|^ 



(E19) 



Hence 



riijrriji af3 

= E (^ m j\ s ±\j\ m r)(f m j'\^ m ^\ a ) 

mjmj/msm'^a^ 

x(3m 3 ; ^a\a*p\3m 3 ; -0) (3m 3 ; ^(3\jm 3 ) 
xY 3 * m3 (r)Y 3m , 3 (r) (E20) 
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where 



(jm^S^jm'j) = Y U m j\ 3m 3^ a )( 3m ^l0\j m 'j) 



x(3m 3 ; ^a\S±\3m' 3 ; -(3) 



for j = 5/2; = — o m >. mjTl 

for j = 7/2; = — o m '. m Tl . 



14 



(E21) 



Then 



= E^I S ±b'^ T 1)0^ T 1|3to 3 ; ~ T ~) 



x(3m 3 :±± hjm^Y^.^if)]' 



1 1. 



1 1, 



= \(jmj T l|3m' 3 ; ^ T ^)| 2 |(3m 3 : ^ ± 



x\Y 3m]Ti (?)\ 2 - 
1) contribution from j = 5/2. Defines 

2 2 2 2 2 2 



2 ~ 2 IJ "^ /I 
(E22) 



14 



x(±) 



7±2(m J Tl) |y i(f)|2 



v 7±2(m jT l)7T2m 



14tt 



l + 2cos 2 6}=g M fl 



i± 
1 1 



Then 



/j>J>(0) = — (1 - 8 sin 2 + 6 sin 4 



2 2 

2)Contribution from j = 7/2. Define 
/£r(0) =/#+/?* 

2 2 2 2 2 2 



(E23) 
(E24) 



14tt 



[5 + 3 cos 2 0] = g 1= if\t (E25) 



Then 



fi 7 (0) = — \2 + 12 sin 2 0-9 sin 4 i 



(E26) 



3) Contribution from Van-Vleck terms. Define 
fii(0)= Y {j,m j \S ± \i'm r \)(j'm j .\f)a^ f} (f\j 1 m,j' / 

2 2 * ' ,H 

rrij rrijf a/3 

Y ti> m j\ s ±\f> m j')ti' m j'\3 m 3>~ a ) 



x(3TO 3 ;ia|cr^ /3 |3m 3 ;i/3)(3m 3 : \f3\jm ) 



2 1 a 'P 1 2 



(E27) 



with 



(j = 2> m i\ s ±\j' = \ m r) 



5 11 1 

E (9 TO 'l 3m 3; -a)(3m 3 ;-a\S±\3m' 3 ;-f3) 



1 7 

x(377z 3 ;-/3|-m) 

T^V(5T2m)(7T2m) 



(E28) 



where to' = m ± 1. 



Insert in, the Eq. ( [E28| ) into the Eq. ( |E27D gives 



,± ^^ V(7 + 2to)(5 + 2to) / 7 ,111. 
/s? = VT 71 (-to|3to± -; - T 



14 

1 1 15 
— : — ± — — 7 

2'2 2 1 2 



2' 2 



x(3TO±^;^±^|>±l)|r 3m± i(r)| 2 



Then 



[2 + sin 2 ( 

147T 



3 sin 4 01 



(E29) 



(E30) 



For cubic crystal, CeSn 3 , we fix the field in the z direction 
and do the reference frame averaging over the cases of the 
Sn atom in the xz plane or the xy plane. Note that yz 
plane is equivalent to the xz plane. 

For the YbCuAl case which has the hexagonal symme- 
try, we have to consider three possible field directions, 
along the x, y and z axis. Then J z — cos aJ z — sinaJ x 
and a z = Aa z — B<r T . where 



A= i(25cos 5 - 26cos 3 + 5cos0) 
4 

B = - sin 0(25 cos 4 9- 14 cos 2 + 1) 
4 



(E31) 



which is derived below. The explicit form of the rotation 
matrices is given by 



E 



(j + to — k)\(j — k — m')\(k — to + m')\k\ 



X (COS l)2j-2k+m-m' ^ _)2fc-m+m' 



(E32) 



where we take the sum over k whenever none of the 
arguments of factorials in the denominator are nega- 
tive. Al has an / = 5/2 nuclear spin and the NMR 
shift was obtained from derivative spectra of the central 
(1/2 <-> -1/2) NMR transition. 

d[^(d)=d ii) 1 _ 1 (6) 

2 2 2 2 
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^ (-l) fc 3!2! . 0. 5 _ 2k . . 0. 2fc 
= E (3 -fc)!(2-*:)!ib!fc! (COS 2 ) ^ 

1 

= - cos - [5 cos 2 — 2 cos — ll = a 

2 2 L J 

d ( 1 f) 1 (0) = -d (| 1 ) iW 

2 2 2 2 

2! 3! 



(E33) 



1 

= --sin-[5cos 2 + 2cos0- 1] = -j3 (E34) 
Then A and £? arc defined as 



A = a 2 - f3 2 = - (25 cos 5 - 26 cos 3 + 5 cos (9) 

B = 2a/? = | sin 0(25 cos 4 0-14 cos 2 + 1) (E35) 



Then 



D 



Sa- («, 9)= A cos affr (0) + ^ sin a[/+, (0) + Jj } , (0)] . 

(E36) 

The explicit values for j = 5/2 and j = 7/2 are given by 



/§ | (a, 9) = — [A cos o(l - 4 sin 2 

2 2 Z07T 

-5sina(3 - 2 sin 2 0)] 
2 

/z z (a, 9) = — [2 A cos a(l + 3 sin 2 
2 2 in 

+Bsma(5 + 3 cos 2 0)]. 



(E37) 



(E38) 



For YbCuAl, we don't include the van Vleck contribution 
because of the large spin orbit splitting as discussed in 
section [V F . 



APPENDIX F: INCLUSION OF COHERENCE 



When multiple scattering is accounting for, the con- 



duction electron Green's function in the Eq. (CI ) become 



dressed by self energy corrections which account for mul- 
tiple scattering off at 4/ sites. As a results, I(e k ,e k r) m 



Eq. C2 is generalized to 



Ijj>(e k ,£k') = -^M kr )h(k'r)^Y^ J — 



d£ A c (e k ,£) 

ILO — £ 



7T VjJ — J 7T W — £ 

dp Aj m (p) f dp' Aj, m ,(p') 



7T i(uj + v) — p J 7T i(u> + V) — p' 

(Fl) 

We can define the conduction electron self energy S c by 



1 



(F2) 



Then 



A c (e k ,0 = -ImG c (k,0 
1 



= -Im 



£-e fc -£ c (0' 



(F3) 



The conduction electron self energy is calculated using 
the average-T matrix (ATA) approximation. 
The v summation is the same as Eq. (C3). Let's do the 
uj sum, 



i i 



i 



P w iUJ ~ ^ iUJ ~ £ ( P ~ iUJ ~ ^ 

_ 2/(Q-2/(p-C) 

(e-eoce + c-p)' 

With Eq.(F4), the Ijj/(e k ,e k /) is written 



(F4) 



I j f(e k ,e k ,)=j 3 (kr)j 3 (k'r) / -^ c (e fc ,£) / -^A c (e fe ,,£') 



7T 



r/£' 



7T 



;; (2/(p)-2JV(C))(2/(Q-2/(p-C)) _ (F5) 



(£-£')(p-p')(£ + C-p) 

In the above equation the £' integration is given by 
df ^(eifc'.O 



ReG c (fcV,£). 



(F6) 



Putting, Eq. (JF6J), Eq. (|C5j) and Eq. (|C6|) together with 
/(p- 0(e _/3C + e"^) = e^ 3 ", %'(e fe) e fc ,) becomes 

IiA^k,^) =h{kr)h{k'r)^- f ^-A c (e k , £)ReG c (fcV, £) 
dC / dp r (B(-) (C) A,- m (p)ReG jW 



IT J IT £ + C — P 

J B(C)4™ ) (p)ReG J w(p)(/(0-l) 1 



c + c-p 



(F7) 



Where J(£) is defined by 



j(0 = / rffc 



k 2 j 3 {kr) 
'£-e fe -E c (0' 



(F8) 



We can write 



fc 2 dfc / k' 2 dk'I jf (e k ,e k -) 



^ImJ(OReJ(0[/(0^'(0 + (1 - /(O)WOL 

7T 

(F9) 



where 
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*vAQ = ^ -=■ 



2 fdC fdp(B(-\OA jm (p)ReG j/m/ (p) 



_L f ^ B (-)( C )ReG. m ( C + OReG J v ro ,(C + 

(F10) 



4/ J 7T ./ 7T 

4/ . 



2 /"dC fdp B(QA { j J(p)BeG j , m '(p)) 

Z4fJ K J 7T 



^- 1 ^^- ) (p)ReG /m ,(p)Rc^(p - 



(Fll) 
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FIG. 2. Schematic configuration diagram for the Uranium 
compound. f 2 , j = 4 and / 3 , j = 9/2 spin-orbit states are 
split by the crystal electric field. All the notations are ex- 
plained in the text and all the values are listed in Table |n|. 



energy 



1 


5 





r 

/' j=7/2 
/' j=5/2 



density of states 

conduction band 



impurity states 



FIG. 1. Schematics of the U = oo single impurity Ander- 
son model for Ce ions. D is the conduction band width and 
V is the single particle hybridization strength. The on site 
Coulomb interaction U is infinity in our model. 
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FIG. 3. Leading order Feynman diagrams for self energies, 
(a) pseudoboson self energy, £0(2)1", (b) pseudo fermion(/ 1 ,j 
multiplet) self energy Y,j m (z). 
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FIG. 4. Leading order Feynman diagram for static mag- 
netic susceptibility. In this diagram, only / electrons are cou- 
pled to the field. 



FIG. 5. Feynman Diagram for the Knight shift calculation 
for Ce ions. This is the lowest order diagram of coupling 
between Ce local moment and nuclear spin in the infinite U 
Anderson model. All the propagators are explained in the 
figure. For the incoherent calculation the conduction electron 
propagator is a bare electron propagator and when the the 
coherence effect for conduction electron is included ( multiple 
scattering off of /-sites) , it becomes a dressed propagator. 
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FIG. 7. Knight shift K(r,T) at fixed angle = and 
Kondo temperature To = 430^ for Ce ions on a fine scale. 
We use the dimensionless variable kpr with the Fermi wave 
vector kF = 0.65A -1 instead of r. Around kpr = 3.2 the lines 
of T = 104K and T = 52K are crossed. This can be explained 
that the Knight shift summation converges faster at higher 
temperatures and it has longer range at lowpi temperatures 
as mentioned by E.S. S0rensen and I. AffleckE5c3. 



10.0 



FIG. 6. Feynman Diagram for the additional contribu- 
tion to— tl»ei Knight shift which was noted by S0rensen and 
Afflec]<EiM This is significant in magnitude for distances 
beyond the Kondo screening cloud radius and not included in 
our calculation. This contribution is of potential importance 
for any dilute syste m, bu t is not important in our lattice cal- 
culation ( see, Sec. [II D| for discussion.). All the propagators 
are explained in the Fig. ^| 



5.0 - 



E ^ 
U 



0.0 



-5.0 1 1 1 1 1 1 

0.0 3.0 6.0 9.0 12.0 15.0 18.0 
KFF 



25 



FIG. 8. Knight shift K(r,T) at fixed angle = 
and Kondo temperature To = 430K for Ce ions. We 
use dimensionless variable Icft with the Fermi wave vector 
fci? = 0.65A -1 instead of r. This shows an oscillatory RKKY 
like behavior. The amplitude is decreased as the temperature 
is increased and the distance r is increased. The calculations 
are done at the same temperatures in Fig ^. 
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FIG. 10. Knight shift K(r,T) at fixed angle 9 = tt/2 
and Kondo temperature To = 430K for Ce ions. We 
use dimensionless variable kpr (with the Fermi wave vector 
kF = 0.65A -1 ). This shows an oscillatory behavior like the 
RKKY interaction form and the amplitude is decreased as the 
temperature is increased and the distance r is increased. The 
calculations are done at the temperatures in Fig [|. 



FIG. 9. Knight shift K(r,T) at fixed angle 9 = tt/2 and 
Kondo temperature To = for Ce ions on a fine scale. We 

use dimensionless variable fc^r (with the Fermi wave vector 
k F = 0.65A- 1 ). The Knight shift has maximum amplitude 
around T = W4K. Around k F r = 2.8 the lines of T = 174ff 
and T — 52K are crossed. 
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FIG. 11. Calculated Knight shift K(T) for a single Ce site 
vs. susceptibility x(T) for varying separation with the Kondo 
scale used to fit the CeSn3 x(T) (see Fig. ^) . For each plot the 
angle is held at 9 = with respect to the nuclear moment-Ce 
axis. The diagram of Fig. | is used to calculate K(T). The 
magnitude of the non-linearity diminishes as is reduced. 
The theoretical Knight shifts have been shifted by offset and 
scale factors to match the susceptibility; this does not affect 
the relative size of the anomaly. 
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FIG. 13. Temperature dependence of Sn Knight 
Shift K(Tfaand Ce xCO (both calculated and experimen- 
tal resultajlj) for CeSn3. The theoretical K(T) is calculated 
using the diagram of Fig.^ and with To chosen to fit the exper- 
imental x{T) data. A full (incoherent) lattice sum is carried 
out over several hundred shells of atoms. 
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FIG. 12. Calculated Knight shift K(T) vs. susceptibility 
X(T) for a single Ce site at fc^r = 3.33 from a nuclear moment 
and angle O = 0. Fixing the /-level energy £5/2 = — 2eV, 
and the spin orbit splitting A ao = 0.29eV, the hybridization 
is varied to illustrate the dependence of the nonlinearity on 
the magnitude of To (which ranges from 750 K to 130 K in 
these calculations as V varies from 0.165el^ to O.lSOeV). The 
diagram of Fig. ^ is used to calculate K (T). The magnitude of 
the non-linearity diminishes as To is reduced. The theoretical 
Knight shift has been shifted by a common offset and scale 
factor to match the susceptibility. 
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FIG. 14. Temperature dependence of Knight Shift K(T) 
with small constant damping — ImE/ = 0.01 in incoherent 
lattice sum. In this figure K(T) = J2k F r<k F r K ( r i T l. The 
calculations are done at the same temperatures in Fig. m The 
small absolute magnitude corresponds to the high tempera- 
ture and it increases as the temperatures goes down. The 
convergence is slower than the coherence lattice sum. 
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FIG. 16. Temperature dependence of Knight Shift 
K(T) with coherent lattice sum. In this figure 
K(T) = y~\ , K(r,T). The calculations are done at 
the same temperatures in Fig. M. The small absolute magni- 
tude corresponds to the high temperature and it increases as 
the temperatures goes down. The convergence is faster than 
the small damping case but slower than larger damping case. 



FIG. 15. Temperature dependence of Knight Shift K(T) 
with large constant damping —ImE/ = 1 in incoherent lattice 
sum. In this figure K(T) = Y,k F r<k F r K ( r > T )' The calcula- 
tions are done at the same temperatures in Fig. 0. The small 
absolute magnitude corresponds to the high temperature and 
it increases as the temperatures goes down. They converge 
very fast, an impurity at large distance doesn't contribute to 
the total Knight shift because of the large damping, i. e. short 
life time of the conduction electrons. 
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FIG. 17. Temperature dependence of 119 Sn Knight Shift 
K(T) with coherenc&ieffect and Ce x{T) (both calculated and 
experimental resultsau) for CeSn3. The theoretical K(T) is 
calculated using the diagram of Fig. |^ and with To chosen to 
fit the experimental xCO data. A full coherent lattice sum 
is carried out over several hundred shells of atoms using the 
average T-matrix approximation to approximate the coherent 
conduction self energy. 



50 
00 
50 
00 
50 
00 



0) 
O 

62 



fi 2 

^1 

n 
O 



t 



a susceptibility x(exp) 
► /isr Knight shift(exp) 
/j,st Knight shif t(incoherent) 
,usr Knight shift(coherent) 



CeSna 



0.50 



00 200 500 400 500 600 
T(K) 



o 



2 
a 

CO 

O 



5.50 
2.50 
1 .50 
0.50 
-0.50 
-1 .50 
-2.50 



A o<x>»< 







CeSn3 



* & susceptibility x( ex p) 

► /.isr Knight shift(exp) 

/^sr Knight shif t(incoherent) - 

* wsr Knight shif t(coherent) 



00 200 500 400 500 600 
T(K) 



FIG. 19. Temperature dependence of positive muon 
Knight Shift K(TXsxLd Ce x( T ) (both calculated and ex- 
perimental resultffla) for CeSn 3 . The theoretical K(T) is 
calculated assuming that muon sits at the middle of the axis 
between the Ce atoms. Other parameters are same as for the 
119 Sn NMR Knight shift calculation. 



FIG. 18. Temperature dependence of positive muon 
Knight Shift K(TAr-e.nd Ce x( T ) (both calculated and ex- 
perimental resultsflO) for CeSn 3 . The theoretical K(T) is 
calculated assuming that muon sits at the center of cubic 
unit cell. Other parameters are same as for the 119 Sn NMR 
Knight shift calculation. 
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FIG. 20. Temperature dependence of positive muon 
Knight Shift K(TAr-e>nd Ce x( T ) (both calculated and ex- 
perimental resultsflo) for CeSn 3 . The theoretical K(T) is 
calculated averaging the ^sr Knight shift from both positions 
with fractional occupancy ratio / = 2/3(2/3of the muon is 
at the center of CeSn3 unit cell.). 
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FIG. 21. Temperature dependence of 27 Al Knight Shift 
A'(T|l|-aad Yb x{T) (both calculated and experimental re- 
sult JitEd) for YbCuAl. The theoretical K(T) is calculated 
using the diagram of Fig. ^ and with To chosen to fit the 
experimental x{T) data. A full (incoherent) lattice sum is 
carried out over twenty four thousand atoms. 
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FIG. 22. Temperature dependence of 89 Y Knight Shift 
A'(T|)-|pjd U x(T) (both calculated and experimental re- 
sultsE30) for Y0.sU0.2Pd3. The theoretical K{T) is calcu- 
lated using the diagram of Fig. |E| and with To chosen to fit 
the experimental x(^) data. A full (incoherent) lattice sum 
and impurity configuration averaging is carried out over a 
several shells. 



TABLE I. Ground states and parameter values for, Ce and 
Yb ions model. j g is the ground state spin- orbit multiplet 
angular momentum and j ex is the excited state multiplet an- 
gular momentum. N g (N,> x ) is the degeneracy of the ground 
(first excited) state multiplet. e/ is the energy of the ground 
configuration, and A ao is the energy difference between the 
j g and jex multiplets 





Ce 


Yb 


Configuration 


4/ 1 


4/ i3 




4/0 


4/ 14 


jg 


5/2 


7/2 


Jex 


7/2 


5/2 


N 3 


6 


8 


N ex 


8 


6 


e f 


-2eV 


-leV 




0.29eV 


1.3eV 



TABLE II. Parameter values for Y0.sU0.2Pd3 in units of 
the conduction electron band width D = 3eV. e/ is the energy 
of j = 4 multiplet and Aij is the energy difference between 
crystal field split Ti and Tj states. F is the single impurity 
hybridization width. For the definition of x and W, see Ap- 
pendix ^. X2 and Wi are for f 2 configuration and X3 and W3 
are for f A configuration. 





-0.333 


A 35 


0.003 


A34 


0.008175 


A31 


0.019621 


Ae8a 


0.000136 


Ae8b 


0.013190 


r 


0.15 


X2 


-0.648 


W 2 


-3.95 x 10" 4 


X3 


0.3693 


W 3 


2.746 x 10~ 4 
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TABLE III. Explicit values of both incoherent and coher- 
ent Knight shift at high temperatures. They have similar 
values because coherent lattice effects are small at high tem- 
peratures. These are unadjusted data. 



Temperature 


K 


night Shift 


Incoherent 


Coherent 


925K 


-0.3579 


-0.44855 


615K 


-2.4078 


-2.42525 


555K 


-3.0157 


-2.84562 


490K 


-3.7184 


-3.34231 


430 


-4.5380 


-3.95714 



TABLE VI. Crystal electric field energy eigenstates for 
j = 5/2 multiplet in the cubic symmetry. 



j = 5/2 Multiplet 


States 


|r 7 5/2J ;T /I) 


-VI l± 5/2) + IT 3/2) 


|rf /2) ;+,!/!} 


1 ± 1/2) 


|r 8 5/2) ;-,T / j) 


v/f |± 5/2) + | T 3/2) 



TABLE VII. Crystal electric field energy eigenstates for 
j = 7/2 multiplet in the cubic symmetry. 



TABLE IV. Crystal electric field energy eigenstates for 
j = 4 multiplet in the cubic symmetry. The coefficient is 
independent of x 



j = 4 multiplet 



states 

~W> yf5\*) + Vt^W + V^I -4) 

|r 3 ;+) 75(|2) + |-2» 

|r s ;-) -v^l 4 ) + ^10)-^ -4) 

|r 4 ;0) _^(| 4 )-|-4)) 

|r 4 ;±i) I +3) + vTl ± 1> 

|r 5 ;0) _^(| 2 )-|-2)) 

|r 6 ;±l) ^| + 3)-V||±l) 



7/2 Multiplet 



States 



|lf /2J ;T /i) 



±7^ | T 7/2)±V^ l±V2) 
l4 7/2) ;T/I) ±^ |±5/2) + i | + 3/2) 



-(7/2). 
8 • 

-(7/2). 



+,T/4) 



; - t/i) 



±1 



+ 7/2) + VtI l± 1/2) 
±5/2) + ^ | + 3/2) 



TABLE VIII. Selection rules for the Anderson hybridiza- 
tion between f 2 F and f V CEF states. The meaning, for 
example, is that a Te conduction electron doublet can com- 
bine with the f 2 Y-j, doublet to make an / 3 Ts quartet. 



TABLE V. Crystal electric field energy eigenstate for 
j = 9/2 multiplet in the cubic symmetry, the coefficient a^'s 
and bj's depend on the x which is the parameter which de- 
pends on the ratio of the fourth and sixth degree cubic field 
in the Hamiltonian of the crystal electric field. In this calcu- 
lation, we use x 3 = 0.3693 and W 3 = 2.746 x 10~ 4 to have 
j = 9/2 Te for the ground state of the / 3 configuration and 
xi = -0.648 and W 2 = -3.95 x 10~ 4 to have j = 4 T 3 for the 
ground state for f 2 configuration. 



j = 9/2 multiplet 



states 



|r 6 ;±) V&liD + VslTD + VisliD 

|r 8a ; 1, ±) ai a | + §) + 02| + 1) + as | ± §)) 

|F 8a ;2,±) 6l |±| ) + &2 | T |) 

|r 8b ;l,±) a 6 |±|) + o 7 | + |) + o 8 |±i)) 

|r 86 ;2,±) b 3 |±f) + &4| + §) 







/ 3 r CEF states 
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r c 


r 6 


»r 3 


r 8 


r 7 1 


»r 3 


r 8 
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r 6 © r 7 © r 8 


r 6 < 


»r 4 


r 6 ffir 8 


r r < 


g>r 4 


r 7 ®r 8 


r 8 - 


g>r 4 


r 6 © r 7 ® 2r 8 


r 6 < 


»r 5 


r 7 ®r 8 


r 7 < 


»r 5 


r 6 ®r 8 


r 8 - 


»r 5 


r 6 © r 7 ® 2r 8 



a For x = 0.3693, oi = -0.1290, a 2 = -0.1582, a 3 = 0.9789, 
bi = 0.7361, b 2 = -0.6769, a 6 = 0.7800, a 7 = -0.6358, 
a s = 0.0016, b :i = 0.6769, and b 4 = 0.7361. 
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